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Abstract. The paper addresses aggregation issues for composite (modular) 
solutions. A systemic view point is suggested for various aggregation problems. 
Several solution structures are considered: sets, set morphologies, trees, etc. 
Mainly, the aggregation approach is targeted to set morphologies. The aggre- 
gation problems are based on basic structures as substructure, superstructure, 
median/consensus, and extended median/consensus. In the last case, prelimi- 
nary structure is built (e.g., substructure, median/consensus) and addition of 
solution elements is considered while taking into account profit of the addi- 
tional elements and total resource constraint. Four aggregation strategies are 
examined: (i) extension strategy (designing a substructure of initial solutions 
as "system kernel" and extension of the substructure by additional elements); 
(ii) compression strategy (designing a superstructure of initial solutions and 
deletion of some its elements); (iii) combined strategy; and (iv) new design 
strategy to build a new solution over an extended domain of solution elements. 
Numerical real- world examples (e.g., telemetry system, communication proto- 
col, student plan, security system, Web-based information system, investment, 
educational courses) illustrate the suggested aggregation approach. 
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1. Introduction 

Traditional approaches to generation of new design solutions are based on mod- 
ifications/improvements of existing products/systems (e.g., [37]). Often one exist- 
ing product /system is used as a basic solution for modification/improvement (e.g., 

nnu, urn nm nn, m, \m, \m, mn) (Fig. 1.1). 



I Existing 
product/ 
V system 

Fig. 1.1. Modification/improvement of existing system 




Modification/ 

improvement 

process 




In this article, the generation of a new modular design solution is examined as 
aggregation of a set of initial existing modular (composite) design solutions. The 
considered aggregation process is depicted in Fig. 1.2: 

Find an aggregated composite solution S a " from a set of initial composite solu- 
tions {Si, S T , S m }, i.e., {Si, S T , —, S m } S a ". 
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Fig. 1.2. Illustration for aggregation process 

Note, the considered aggregation process has to be based on an analysis of engi- 
neering and/or management application(s). An expert-based scheme (framework) 
for problem solving (Fig. 1.3) is used as a basic framework for aggregation of 
composite (modular) solutions. 

Here it is necessary to point out main properties of the expert-based approach: 

1. Each operation/step of the solving process has to be executed while taking into 
account an applied expert-based analysis of the problem situation or subsituation. 

2. At each step of the solving process an expert opinion is the most important 
and the domain expert can correct the solving scheme and intermediate and/or 
resultant solutions. 
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Fig 1.3. Expert-based approach to problem solving 

Generally, two aggregation cases have to be considered: 
Case 1. Usage of elements from m initial solutions, 

Case 2. Usage of elements from m initial solutions and additional design elements 
from an extended design domain as well. 

In case 1 three solving strategies can be examined. First, the evident strat- 
egy consists in an analysis of the set of m initial solutions and building a system 
substructure ("system kernel" K, e.g., a "basic system part" as substructure, me- 
dian/consensus of the initial solutions) that can be extended or modified. Several 
methods can be used to build the "system kernel" K: (a) substructure of the initial 
solutions, (b) median structure (or consensus/ agreement) for the initial solutions, 
(c) extended substructure, (d) extended median structure (or extended consensus), 
(e) selection of the "best" system element for each specified system part. 
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Here expert judgment of expert(s) domain can be used at each stage of the 
solving process. 

Second, a superstructure for set of m initial solutions has to be designed fl 
(e.g., combining the initial solutions, "covering" the initial solutions). Then, the 
superstructure is compressed via deletion of the less important elements. 

In case 2, it is possible to consider the following: 

to build a generalized domain of design alternatives (i.e., solution elements) and 
to design a new solution (including usage of new design elements). 

Here the extended design domain may be obtained via the following two ways: (a) 
covering/approximaiton of m initial solutions by a new design alternative domain, 
(b) examination of a new design alternative domain. 

In the article, the following structures will be examined as basic ones: (i) sets 
(e.g., [S3), (ii) rankings (e.g., [29], [35], [36], [101] ) . (hi) set morphologies (e.g., [9], 

m, mm, urn nnsi, una, cm w tr ees ( e . g ., m, m)- 

General structure (i.e., system model) (A) is: A= ( T,P,D,R,I ), where the fol- 
lowing parts are considered ( |101j . |106) . [108] ): (i) tree-like system model T, (ii) 
set of leaf nodes as basic system parts/components P, (hi) sets of DAs for each leaf 
node D, (iv) DAs rankings (i.e., ordinal priorities) R, and (v) compatibility esti- 
mates between DAs I. Fig. 1.4 illustrates the generalized architecture (structure) of 
examined modular systems/solutions (A). Note the following significant structures 
are considered as well: (a) morphology: $ = ( P.D ); (b) morphological set: $ = ( 
P,D,R,I ); (c) morphological structure (tree): ( T,P,D,R ); (d) morphological 
structure with compatibility (i.e., with compatibility of DAs): <fr= ( T,P,D,R,I ). 




System parts/ 
components 
(leaf nodes P) 



Sets of design Compatibility 
alternatives among 
(morphology D) design 

/ alternatives (I) 




Fig 1.4. Architecture of modular system (A) 

In Table 1.1, building problems for basic substructures/superstructures are briefly 
pointed out. 

Finally, two basis system problems are faced: (a) revelation/ design of "system 
kernel", (b) building an extended design alternative domain. Four solving aggre- 
gated strategies are considered: 

(1) extension strategy: designing "system kernel" based on set of initial solu- 
tions (e.g., substructure, median/consensus) and extension of "system kernel" by 
additional elements; 

(2) compression strategy: designing a superstructure of initial solutions and 
deletion of some its elements; and 
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(3) combined strategy (extension, deletion, and replacement operations for sys- 
tem elements over a preliminary aggregated solution), and 

(4) design strategy: building an extended design domain of solution elements 
and designing a new solution. 

The above-mentioned strategies are based on combinatorial models (as underlay- 
ing problems): multicriteria ranking/selection, knapsack problem, multiple choice 
knapsack problem, combinatorial morphological synthesis. 

The suggested approaches are illustrated through numerical examples including 
applied examples (e.g., telemetry system, communication protocol, hierarchical se- 
curity system, plan of students art activity, combinatorial investment, Web-based 
information system, notebook, educational courses). 



Table 1.1. Basic structures and bibliography references 



Initial structures 


Target structure 


References 


Sets 


Subset, superset 


0, [77] 


Strings/sequences 


Common subsequence, 
common supersequence, 
median string, consensus 


®, 0, eh. m m, m, m, 

[HE [S3, G3I, El, [B2], [Sj, 

[126], [129], [130], [135], |143|, 
[145], |151j. [153] 


Rankings 


(i) Consensus/median 

(ii) Fuzzy ranking 


ELI, HE Ei], [33, [S3, nag, 

[72], [88], [89] 

[71], [91], [100], |1()1| 


Trees 


(i) Agreement / consensus 
tree 

(ii) Agreement forest 
(hi) Supertree 


0, m, M, [S2], [54], [76], 
[133], [146] 

|30j. |66|. |138]. [164] . [165] 
PI, P|, [75], [142] 


Graphs 


Common subgraph, 
common supergraph 


[53]. |81j. |104j. [106] 



2. Auxiliary Problems and Aggregation Strategies 

2.1. Basic Auxiliary Problems. In the case of two initial solutions as element 
sets Ai and A2, Fig. 2.1 illustrates substructure Sa 1 a 2 Q (^l&^b) and superstruc- 
ture SU1A2 ^ (AiUA 2 ) (via Venn-diagram). 



Superstructure Sa 1 a 2 — (^iU-^2/) 

Substructure <SUiA 2 ^= {A1&A2) 




Fig. 2.1. Illustration for substructure, superstructure 
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Further, let S = {Si, Si, S n } be a set of initial solutions (structures). Let 
function p(Si 1 , Si 2 ), i\,i2 £ {1, 2, n} be a proximity or a metric for the solutions. 
The following main basic auxiliary problems can be examined: 

Problem 1. Find a maximum substructure: 

n 

S = arg max (\S'\), VS" e f] Si. 
{s ' } i=i 
Problem 2. Find a minimum superstructure 

71 

S = arg min (|S*"|), VS*" e (J 
{S " } i=i 
Now let us consider definitions of medians for the above-mentioned set of initial 
sets S (e.g., @U], [HI], USD], Q35]): 

(a) median ("generalized median") M 9 is: 

n 

M 9 = arg miriMeD p(M, A)), 

i=l 

where D (£> 3 S) is a set of structures of a specified kind (searching for the median 
is usually NP complete problem); 

(b) simplified case of median (an approximation) as "set median" M s over set 

S: 

n 

M s = arg miriMes (^ p( M ^i)). 

i=l 

Here a representative from S = {Si, Si, S n } is searched for. Computation of 
proximity p(M, Si) is usually NP-complete problem as well. Note a similar "closest 
string problem" is widely applied in bioinformatics (e.g., [S5], [55], [53], [57], |122j . 
[125] . |163j V Finally, the following problem 3 and problem 4 can be considered: 

Problem 3. Find "set median" M s . 

Problem 4- Find "median" ( "generalized median" ) M 9 . 

Problem 5. Find an extended median/consensus structure via addition to (or 
correstion/editing of) the basic median/structure some elements while some re- 
source constraint (s). Here some elements are added to the median set (problem 3 
or 4) while taking into account profit and required resource for the addition. 

The problems are considered for the following kinds of structures: (i) sets, (ii) 
set morphologies, (iii) trees, and (iv) trees with set morphologies. In the case of 
vector- like metric/proximity p{Si 1 , Si 2 ), Pareto-efficient solutions are searched for 
in problem 3, problem 4, and problem 5. 

2.2. Building of "System Kernel". The basic auxiliary problem consists in de- 
signing the "system kernel". Several methods can be used to build the "system 
kernel" K: 

(i) substructure of the initial solutions, 

(ii) median structure (or consensus/agreement) for the initial solutions, 



AGGREGATION OF COMPOSITE SOLUTIONS: STRATEGIES, MODELS, EXAMPLES 7 

(iii) extended substructure (or extended median structure, extended consensus), 
and 

(iv) a two-stage framework: (a) specifying a set of basic system parts (as a subset 
of the system parts/components), (b) selection of the "best" system elements for 
each specified basic system part above. 

Sometimes, a subsolution (i.e., "system kernel" as substructure) is a very small 
subset. In this case, it is reasonable to use a special (more "soft") method to select 
elements for "system kernel". Let S = {Si, Si, S n } be a set of initial solutions 
(structures). Then, element e will be included into (added to) "system kernel" if 
r]i > a (e.g., a > 0.5) where r\ e is the number of initial solutions which involve 
element e. The usage of this rule will lead to an extension of the basic method for 
building the system substructure. 



2.3. Aggregation Strategies. Four basic aggregation strategies are examined: (i) 
extension strategy (extension of "system kernel"), (ii) deletion strategy (compres- 
sion of a superstructure for initial solutions), (iii) combined strategy (i.e., extension, 
deletion, and replacement operations for system elements over a preliminary aggre- 
gated solution), and (iv) design strategy (new system design). 
The extension strategy is the following: 

Type I. Extension strategy: 

Phase 1.1. Analysis of applied problem, initial solutions, resources, solution 
elements. 

Phase 1.2. Revealing from m initial solutions a basis as a subsolution or "system 
kernel" K (e.g., subset of elements, substructure, median/consensus). 

Phase 1.3. Forming a set of additional solution elements which were not included 
into the basis above. 

Phase 1.4- Selection of the most important elements from the set of additional 
elements while taking into account the following: (i) profit of the selected elements, 
(ii) total resource constraint (s), (iii) compatibility among the selected elements and 
elements of the basis (i.e., "system kernel"). 

Phase 1.5. Analysis of the obtained aggregated solution(s). 

Fig. 2.2 illustrates the extension strategy. 




Fig 2.2. Scheme of extension strategy 
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The second (compression) strategy is based on preliminary union of all elements 
of m initial solutions and deletion of the non-important elements to satisfy some 
resource constraint(s) (Fig. 2.3): 

Type II. Compression (deletion) strategy: 

Phase 2.1. Analysis of applied problem, initial solutions, resources, solution 
elements. 

Phase 2.2. Union of all elements from m initial solutions to form a basic super- 
solution (e.g., superset of elements, superstructure). 

Phase 2.3. Forming a set of element candidates to delete. 

Phase 2.4- Selection of the most non-important elements from the set of element 
candidates while taking into account the following: (i) integrated profit of the 
compressed solution, (ii) total resource constraint (s), (iii) compatibility among the 
selected elements of the compressed solution. 

Phase 2.5. Analysis of the obtained aggregated solution(s). 



I 



Initial solutions {Si, .., S T , S m } 

I 



^'SysteiriX /Supersolu 
V kernel" K J Lion ft 



Element-candidates 
for deletion A C Q\K 



I 



Compression strategy (type II) 
1 

Aggregated solution S a " 



Fig 2.3. Scheme of compression strategy 

The third strategy is consists in possible combination of addition operations, 
deletion operation, and correction (replacement) operations. The third strategy is 
(Fig. 2.4): 



Initial solutions {Si, .., S T , S m } 




f 'Systcrm ( Super-\ 
kernel" sotion 


Candidates 
for addition 
W 


Candidates 
for deletion 
B 


Candidates 
for replacement 
C = \(a u ,b v )} 




I 




I 



Extension/Compression strategy (type III) 
1 

Aggregated solution S a " 



Fig 2.4. Scheme of extension/compression strategy 
Type III. Combined strategy: 
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Phase 3.1. Analysis of applied problem, initial solutions, resources, solution 
elements. 

Phase 3.2. Revealing from m initial solutions a basis as a subsolution or "system 
kernel" K (e.g., subset of elements, substructure, median/consensus). 
Phase 3.3. Forming the following: 

(a) a set of additional solution elements as candidates for addition W which were 
not included into the basis above (|WP)if| = 0), 

(b) a subset B C K as candidates for deletion, and 

(c) a set of element pair C — {(a u , b v )\(a u € K)k(b v GK)}, thus a set of correc- 
tion operations is considered as replacement of a u by b v . 

Phase 3.4- Selection of the most important operations (including element ad- 
dition, clement deletion, and element replacement) while taking into account the 
following: (i) profit of the operations, (ii) total resource constraint (s), (hi) compat- 
ibility among the selected elements in the resultant solution. 

Phase 3.5. Analysis of the obtained aggregated solution(s). 

The fourth strategy is targeted to usage of additional elements from an extended 
design domain ( "design space" ) . Here the resultant aggregated solution may involve 
elements which were not belonging to m initial solutions. The fourth strategy is: 

Type IV. Strategy of extended design domain: 

Phase 4.1. Analysis of applied problem, initial solutions, resources, solution 
elements. 

Phase 4.2. Extension of the union of all elements from m initial solutions to 
form an extended design element domain ( "design space" ) . 

Phase 4-3. New design of the composite solution over the obtained design el- 
ement domain while taking into account the following: (i) profit of the designed 
solution, (ii) total resource constraint (s), (hi) compatibility among the selected (i.e., 
resultant) elements of the compressed solution. 

Phase 4-4- Analysis of the obtained aggregated solution(s). 

Fig. 2.5 illustrates the fourth strategy. Given three initial solutions Si, S2, and 
S3. The resultant aggregated solution S a " can involve elements from solution Si, 
S2 and elements from the extended design domain. 



Fig. 2.5. Illustration for aggregation design strategy 

3. Examined Structures, Substructure, Superstructure 

Generally, the following basic kinds of sets are examined: sets, multisets, lists, 
complete lists, and trees (e.g., [43], [93], [147] ). In this material, our main examina- 
tion is targeted to special kinds of composite structures: morphological structures. 
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3.1. Sets. Evidently, sets are basic structures. Fig. 3.1 depicts a numerical exam- 
ple of two initial sets A and B. 



Set A 



• • • • • 

1 3 4 5 6 



• • • 

101112 



Set B 



1 2 



6 7 



9 10 



1*2 



Fig. 3.1. Illustration: two initial sets 
Fig. 3.2 contains examples of subset Sab and superset Sab- 

Subset Sab Superset Sab 



4 



10 



1*2 



1 2 3 4 5 6 7 8 9 101112 



Fig. 3.2. Illustration: subset and superset 
Now let us consider the case of m sets. Let {A\, A+, A m } be initial sets. 
Then, a subset is: S^aj Q DfcLi ^i- A superset is: SiA t } 2 U2=i ^i- 

3.2. Rankings (Layered Sets). Here ranking is examined as a layered set. Let A = 
{1, i, n} be a set of elements/ items. Ranking (a partial order/partition of set 
A) is considered as linear ordered subsets of A (Fig. 3.3): A — UfcLi l-A(fci) (~1 

A{k 2 )\ = if fci ^ k 2 , i 2 ^ ii Vii G A(fei), Vi 2 e ^(fca), fei < fc 2 . Set A{k) is 
called layer fc, and each item i G A gets priority r, that equals the number of the 
corresponding layer. The described partition of A is be called as partial ranking, 
stratification, layered set (e.g., [17], [42], [98], [l00], QH], [106], [O], [139], [T73] V 
Evidently, a linear order of elements from A is ranking as well. Many years ranking 
problems (or sorting) have been intensively used and studied in various domains 
(e.g., PU, [23], [42], [29], [88], [89], [106], [IH], [139], [173], [III). 



Initial elements 
A = {1, ...,n} 



-4i 
TTTT 



I I I I 

TTTT 
I LLi 

Fig. 3.3. Scheme of ranking 

3.3. Multisets and Morphological Sets. Now let us consider two kinds of basic 
structures: multiset and morphological set. 

Let U be a universe (collection of all relevant items) and N be a set of non- 
negative integers (e.g., N = {1, ...,t, ...,m} ). Formally, multiset is a mapping: 



T : U 



N. 



Clearly, three basic operations can be considered for sets and multisets: 
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1. unions: A[j B (for sets A and B, analogically for multisets), 

2. intersections: Af]B (for sets A and B, analogically for multisets), and 

3. complements: A\B (for sets A and B, analogically for multisets). 

A macroset is a (finite or infinite) set of multisets over a finite alphabet. 
Generally, morphological set is defined as follows: 

Definition 3.1. Morphological set is a structure consisting of: 

1. a finite set of integers N = {1, r, m} (each integer r corresponds to a 
system part); 

2. set of elements (alternatives) for each system part r: 

A T = {A T \, A T £, A T q T }, where A T £ is a design alternative; 

3. preference relation over elements of A T (or estimates upon a set of specified 
criteria or resultant ordinal priorities for each alternative p(A T ^)); and 

4- weighted (by ordinal scale) binary relation of compatibility for each pair of 
system parts (a,/?) <E N over elements of alternative sets A a ,Ap: Rj\. a ,Ag- 

A composite system consisting of m parts is examined (Fig. 3.4) ( |101] . [IDB] . |110| V 



System 



Part / 


Part 






r 




Ui) 


■■ a 


Set of 


Set of 




alternatives 
for part 1 



alternatives 
for part r 



Set of 
alternatives 
for part m 



Fig. 3.4. Illustration for multi-part system 



Thus, system morphology Sa (or morphology $ ) is defined as follows A 
{A 1 ,...,A T ,...,A m } (Fig. 3.5). 



System morphology A 




Fig. 3.5. Illustration: system morphology 



Further, numerical examples for two morphologies are presented: morphology 
A 1 (Fig. 3.6), morphology A 2 (Fig. 3.7), sub-morphology Sa 1 ^ 2 (Fig- 3.8), and 
super-morphology Sai,a 2 (Fig- 3.9). 
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Morphology A 1 





■^Ay 








Al2 






A13 




A 23 




Morphology A 2 




Fig. 3.6. Example of morphology A 1 
Sub-morphology S A i A 2 



Fig. 3.7. Example of morphology A 2 
Super-morphology S^.a 2 





Fig. 3.8. Sub-morphology Sa^,a 2 



Fig. 3.9. Super-morphology Sa 1 a 2 



In addition, tables of ordinal compatibility estimates have to be considered. The 
aggregation of estimates (by element) can be based on a special operation (e.g., 
minimal value of element). Table 3.1 presents an illustrative numerical example for 
compatibility estimates and their aggregation. 

Table 3.1. Aggregation of compatibility estimates for A 1 , A 2 : nodes A\,A 2 

Super-morphology Sa u a 2 



Morphology A 1 


Morphology A 2 




A 2X 


^22 


^23 


A 24 




A 2 i A 23 A 2 4 




A 22 A 23 


An 









2 


An 


1 2 


A 12 


2 3 


A 12 


3 


2 


CD 


1 


Al2 


3 2 1 


A 13 


1 2 


A 13 


2 


1 


(D 


1 


Ai 3 


2 3 1 


Al4 


2 2 


A 2i 




2 


2 





3.4. Trees and Morphological Structures. Now let us consider trees and mor- 
phological structures (as morphological trees). Fig. 3.10 illustrates a numerical 
example for two trees T\ and T2. 




Fig. 3.10. Illustration for two initial trees 



Fig. 3.11 depicts corresponding examples of supertree T and subtree T. 
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Supertree T Subtree f 




Fig. 3.11. Illustration for supertree and subtree 

Now let us consider trees with morphologies. Fig. 3.12 illustrates two trees with 
morphologies 9i and 2 . 




Fig. 3.12. Illustration for two initial trees & morphologies 



Fig. 3.13 illustrates examples of supertree and subtree with morphologies 6 and 

e. 



Supertree & morphology 6 Subtree & morphology 9 




Fig. 3.13. Illustration for supertree and subtree & morphologies 



Evidently, tables of compatibility estimates and their aggregation can be consid- 
ered here as well. 

4. Preliminary Illustrative Example for Notebook 

The preliminary illustrative applied example is targeted to representation and 
aggregation of three initial personal computers (notebooks) . The considered general 
morphological structure of the notebook is presented in Fig. 4.1: 

0. Notebook S. 

1. Hardware H: 

1.1. Basic computation C: 

1.1.1. Mother board B: B x (P67A - C43(B3) ATX Intel), B 2 (MSI 870A-G54 
ATX AMD), B 3 (ASRoot P67 EXTREME 4(B3) ATX Intel); 
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1.1.2. CPU U: Ui (Intel Pentium dual-core processor T 2330), U 2 (Celeron 
dual-core processor 2330), Us (Intel core 2 T 7200); 

1.1.3. RAM R: Ri (1 GB DDR A-DATA), R 2 (2 GB DDR2 KINGSTON), R 3 
(2 GB DDR3 A-DATA), R 4 (2 GB DDR3 HYPER X KINGSTON); 

1.2. Hard drive V: V 1 (100 GB HDD), V 2 (120 GB HDD), V 3 (160 GB HDD), 
Vi (200 GB HDD); 

1.3. Video/graphic cards J: J x (NVIDIA GeForce CTS 300M), J 2 (GT 400M 
Series), J 3 (ATI Radion HD 5000 M Series); 

1.4. Communication equipment (modems) E: E\ (Internal Modem & Antenna), 
E 2 (None). 

2. Software W: 

2.1. Operation system and safety Y: 

2.1.1. OS O: Oi (Windows XP), 2 (Windows Vista); O3 (Linux). 

2.1.2. Safety software F: Fi (Norton Antivirus), F 2 (Antivirus Kaspersky). 

2.2. Information processing and Internet /: 

2.2.1. Data support and processing D: D\ (Microsoft Office), D 2 (None). 

2.2.2. Internet access (browser) A: A\ (Microsoft Internet Explorer), A 2 (Mozilla); 
A3 (Microsoft Internet Explorer & Mozilla). 

2.3. Professional software Z: 

2.3.1. Information processing (e.g., engineering software) G: G\ (Matlab), G 2 
(Lab View); G 3 (MatCad), G 4 (None); 

2.3.2. Special software development environment P: P\ (C++), P 2 (JAVA); P 3 
(Delphi), P 4 (None); 

2.3.3. Special editors L: L\ (LaTcx), L 2 (None); 

2.3.4. Games Q: Q x (Tctris), Q 2 (Solitaire), Q 3 (Chess). 




Fig. 4.1. Simplified structure of notebook 



Now a simplified illustrative example for four initial solutions Si, S 2 , S3, and ^4 
is examined (Fig. 4.2, Fig. 4.3, Fig. 4.4, Fig. 4.5). Here the tree-like structure is 
not changed and only leaf nodes are considered. Compatibility estimates between 
design alternatives are not considered. A framework of aggregation process for four 
notebooks is depicted in Fig. 4.6 (including two alternative methods to build the 
"system kernel"). 
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M000000000000 

Fig. 4.2. Structure of notebook Si 

M000000000000 

Fig. 4.3. Structure of notebook S2 
B U R V J E O F D A G P L Q 

Fig. 4.4. Structure of notebook S3 



B U R V J E O F D A G P L Q 

Fig. 4.5. Structure of notebook S4 

Substructure for the considered solutions S is depicted in Fig. 4.7, superstructure 
S is depicted in Fig. 4.8. Further, let us consider "system kernel" K as an extension 
of substructure S (Fig. 4.9). 



(^Notebook sT) (^Notebook S3) 



Notebook Si 



Notebook S4 



Building "system kernel" K (via subso- 
lution, via selection of the best system) 
element for each system components) 



(^"System kernel" 



Correction of "system kernel" K (e.g., 
element addition, element replacement) 



I 



Aggregated solution S a " 



Fig. 4.6. Framework of aggregation process 
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Fig. 4.7. Substructure of solutions S 



B 

ft 

B 2 



U R V J 



E 



O F D A G 



U 2 
U 3 



R2 



J2 



O2 



F 2 



ft 

A 3 



Gi 

G4 



ft 

P 3 
Pa 



Q 



Fig. 4.8. Superstructure of solutions S 
B U R V J E F D A G P L Q 
V< II II /•' II I II h ||/ II" II II /' II h II " || * 11^: 






Fig. 4.9. "System kernel" K 

Further, the aggregation strategy as modification of "system kernel" K can be 
applied. A set of candidate modification operations are the following: 

1. addition operations: 1.1. addition for U: U\ or U 2 or U 3 , 1.2. addition for 
F: F x or F 2 , 1.3. addition for P: P 2 or P 3 or P 4 ; 

2. correction operations: 2.1. replacement B\ =>■ B 3 , 2.2. replacement V3 => V4, 
2.3. replacement A\ ^> A 3 . 

Evidently, it is reasonable to evaluate the above-mentioned modification opera- 
tions (e.g., cost, profit) and to consider an optimization model. Later, correspond- 
ing optimization problems (e.g., knapsack problem, multiple choice problem) will 
be examined. An example of the resultant solution S 1 (modification of "system 
kernel" K) is shown in Fig. 4.10. 




Fig. 4.10. Solution S 1 based on modification of "system kernel" K 

On the the hand, building the "system kernel" can be based on multicriteria 
selection and/or expert judgment. Let us consider the following basic structure of 
"system kernel": B, U, R, V, O, F, D, G. For each system component above, it 
is possible to consider a selection procedure to choose the "best" system element 
(while taking into account elements of the initial solution or additional elements as 
well) . For example, we can obtain the following "system kernel" : 

K* = B l -kU 2 -kR 2 *V 3 *E 2 *0 2 *F 2 -kD 2 -kG 2 . 

Further, the system correction process is based on the following operations: 

1. addition: A\, P\, L\\ 

2. deletion: E 2 ; 

3. replacement: Si =>• B 3 , U 2 => Ui, 2 => 0\. 
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A resultant solution S 2 is presented in Fig. 4.11. 




Fig. 4.11. Solution S 2 

5. Metrics and Proximities 

Let us consider similarity measure between objects (in our case: sets, rankings, 
trees, graphs) A± and A 2 . It is often desired that the distance measure (function) 
d(Ai,A 2 ) fulfills the following properties of a metric: 

1. d(Ai,A2) > (nonnegativity) , 

2. = (identity), 

3. d(A u A 2 ) = A = B (uniqueness), 

4. d(A 1 ,A 2 ) = d{A 2 ,A x ) (symmetry), 

5. d(Ai,A 2 ) + d(A 2 ,A 3 ) > d(Ai,A 3 ) (triangle inequality). 

If the function satisfies d(Ai, A 2 ) < 1 it is said to be a normalized metric. 

In many applied domains, the above-mentioned conditions are too restrictive and 
a more weak set of properties is used. As a result, d(Ai,A 2 ) corresponds to more 
weak situations, for example: (i) quasi-metrics, (ii) proximities (without property 
5, e.g., proximity for rankings in [101J. 

Metrics/proximities play the basic role in many important problems over struc- 
tures, for example: approximation, modification, aggregation. There exist three 
basic approaches to similarity /proximity of objects/structures: 

(1) traditional metrics/distances (e.g., [22], [47], [60], [61], [128] ): 

(2) minimum cost transformation of an object/strcuture into another one (edit 
distance) (e.g., [20], [25], [68], [97], [127], [HH], [EH], [m]); and 

(3) maximum common substructure or maximum agreement substructure (e.g., 

®, ®, m, m, \m, cm cm »■ 

5.1. Metric/Proximity for Sets. Let A = {1, i, n} be a set of elements. 
Let us consider ... for two subsets A 1 C A and A 2 C A. The most simple case 
of metric by elements (i.e., distance) is the following: (while taking into account 
assumption {A 1 [j A 2 \ ^0): 

x 2 _ i^pi^ 2 ! 

Pe (A ,A )_!___. 

Further, let Wi G (0, 1] be a weight of element i € A. Then proximity (i.e., metric, 
distance) by element weights is as follows (while taking into account assumption 

EieOl 1 u a 2 ) w i 7^ 0) : 

Pw {A\A 2 ) = l- ^ A1 ^\ W \ 

Now let us consider a simple numerical example (Table 5.1) that involves ini- 
tial set A = {1,2,3,4,5,6,7,8} and two subsets: A 1 = {1,2,4,5} and A 2 = 
{1,2,3,4,5,6,7}. Elements weights are the following: 0.5, 0.6, 0.4, 1.0, 0.7, 0.2, 
0.1, and 1.0. The resultant proximities are: p e (A 1 ,A 2 ) = 1 — | = 0.571 and 
Pw (A\A 2 ) = l-^=0A96. 



18 



MARK SH. LEVIN 



Table 5.1. Illustrative numerical example 



i 


Wi 


A 


A 1 


A 2 


A 1 \JA 2 


a 1 n a 2 


1 


0.5 


★ 


★ 


■k 


■k 




2 


0.6 


-k 


-k 


-k 


-k 




3 


0.4 


-k 




-k 


* 




4 


1.0 


* 


-k 




■k 




5 


0.7 




■k 




-k 




6 


0.2 








■k 




7 


0.1 








-k 




8 


1.0 













In the case of vector weights uii = (w-, wf, w^), a simplified vector distance 
(or proximity) is: 

MA 1 , A 2 ) = (pUA\A 2 ), pL(A\A 2 ), pUA\A 2 )) = 

l^ie(A l UA 2 > w i Lie^LU 2 ^ l^ieiA 1 (j a 2 ) w i 

5.2. Proximity for Strings/Sequences. This problem of proximity analysis is 
important in code design, genom studies, information processing, and mathematical 
linguistics. Generally, three kinds of proximity for strings/ sequences are examined: 
(1) common part of initial strings as substring, superstring (e.g., [7], jB], 00], [55], 
[71], [79], P33]), (2) median string (e.g., [94], Q30]), and (3) editing distance (a 
lenght/cost of an transformation/editing) (e.g., [6], [68], [97], [137j . [160] ). 

Here basic research issues are targeted to the following: (a) complexity analysis 
(e.g., |153j ): (b) design of polynomial algorithms (e.g., [68], [79]; and (c) design of 
approximation algorithms (e.g., [78]). 

5.3. Proximity for Rankings. Here several types of metrics/proximities can be 
used, for example: 1. Kendall- Tau distance [90]; 2. distances for partial rankings 
(e.g., Q3], [H]); 3. vector-like proximity ([98], [100], [TOl]). Further, Kendall-Tau 
distance and vector-like proximity are briefly described (the description is based on 
material from [101] ). 

Let \\gij |[, £ A) be an adjacency matrix for graph G: 

( 1, if*>-j\ 
9ij = I 0, if i — j, 
{ -1, ifi-<j. 

Kendall-Tau distance (metric) for graphs G 1 and G 2 is the following: 

PK{G\G 2 )=Y, \9l 1 -9l I, 

i<j 

where gjj, g 2 ^ are elements of adjacency matrices of graphs G 1 and G 2 accordingly. 

Basic definitions for vector-like proximity are the following ([98], |100) . |101j ). 
Let ^>(S) be a set of all layered structures on A. 
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Definition 5.1. Let us call the first order error Vi G A, and the second order 
error G {A* A\i ^ j VS,Q G as follows: 

5f(S,Q)=7T i (S)-7T i (Q), 

S7.(S, Q) = MS) - irj(S) - (7Ti(Q) - ttj(Q)), 

where ni(S) = / Vi G A(Z) in 5. Thus, for an estimate of a discordance between the 
structures S,Q G ^(S) with respect to i and an integer- valued scale with the 

following ranges is obtained: 

(i) -(to - 1) < r < m - 1 for <5f (S, Q), and 

(ii) -2(m - 1) < r < 2(m - 1) for 5^.(5, Q). 

Definition 5.2. Let 

ar(S ) Q) = ( a5 -( m - 1 ),..., a: - 1 ,a : 1 ,..., a: " , - 1 ) ) 

2/ (5,g) = ( y - 2 (™- 1 ),..., y - 1 , y 1 ,..., y 2 (™- 1 )), 

be vectors of an error (proximity) VS, Q G ^(S) with respect to components i 
(1st order), and the pairs (i, j) (2nd order). The components of above-mentioned 
vectors are defined as follows: 

x r = \{i£A\6?(S,Q)=r}\/n, 

V r = 2\{(i,j) G {A * A\i + JMjiS, Q) = r}\/(n(n - 1)). 

Now denote a set of arguments for the components of vectors x and y as follows: 
O = { — k, negative values as Q _ , and positive ones as f2 + . In addition, we 

will use the vectors x with aggregate components of the following type (similarly, 
for y): 

&2 — k m+1 

x ki,k 3 = x r ,x^~ k = x r ,x- k = J2 x ^ k>0 ' xH =x r + x- r . 

r=fei r=— (m— 1) r=k 

Definition 5.3. Let ^(S 1 , Q)\ = ^2 reQ x r , \y(S, Q)\ = Y, r enV r bc modules of 
the vectors. 

Vector x will be used as a basic one. 

Definition 5.4. We will call vectors truncated ones if 

(1) the part of terminal components is rejected, e.g., 

X ^ S ■ C^) ^ — (^X j X ^ j . . . ^ 3/ ^ 3/ ^ . . . ^ ■ ) • 

and one or both of following conditions are satisfied: k\ < to — 1, &2 < m — 1; 

(2) aggregate components are used as follows: 

x(S,Q) = {x^ kl ,...,x ka -\x ka ' kb ,x kb+1 ,...,x^ k2 ), 

x(S,Q) = (xW,...,x^,...,xW). 

Definition 5.5. Let us call vector x (y): 

(a) the two-sided one, if \ ^ and \ \ = 0; 

(b) the one-sided one, if \Q+\ = or |f2~| = 0; 

(c) the symmetrical one, if — r G 0~ exists Vr G fi + , and vice versa; 

(d) the modular one, if it is defined with respect of definition 5.4 (5.3). 
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Moreover, a pair of linear orders on the components of vectors x and y (definition 
5.2) is obtained: component 1(— 1) -< ... -< component k(—k). 

Definition 5.6. x^S, Q) h x 2 (S,Q), 9,( Xl ) = Q(x 2 ), VS,Q e #(S), if any 
decreasing of weak components x\ in the comparison with x 2 is compensated by 
corresponding increasing of it's 'strong' components (r,p e fl + or — r, -p£ Q~): 



c£ > 0,Vu e 0+(V-u e Q",-r < -u). 



Now let us consider an illustrative numerical example: 

(a) initial set A = {1, 2, 3, 4, 5, 6, 7, 8, 9}, 

(b) ranking S 1 = {A\ = {2, 4}, A\ = {9}, = {1, 3, 7}, A\ = {5, 6, 8}} and 

(c) ranking S 2 = {A\ = {7, 9}, A% = {1,3}, A% = {2, 5, 8}, A\ = {4, 6}}. 
Corresponding adjacency matrices are as follows: 



1^(^)1 - 





-1 





-1 


1 




1 








-1 




-1 


1—1 





1 




-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 





-1 





-1 


-1 


-1 


-1 


-1 


1 


-1 


1 


-1 



-1 1 1 



1 

1 1 

1 



1 

-1 
-1 



-1 . 
1 1 



-1 \ 

1 

-1 
1 

-1 
-1 
-1 
-1 

• J 



\9ij(S 2 )\ = 



( 




i 





1 


1 


1 


-1 


1 


-1 \ 




-1 




-1 


1 





1 


-1 





-1 







i 




1 


1 


1 


-1 


1 


-1 




-i 


-i 


-1 




-1 





-1 


-1 


-1 




-i 





-1 


1 




1 


-1 





-1 




-i 


-1 


-1 





-1 




-1 


-1 


-1 




i 


1 


1 


1 


1 


1 




1 







-i 





-1 


1 





1 


-1 




-1 


V 


i 


1 


1 


1 


1 


1 





1 


• / 



Now it is very easy to calculate Kendall- Tau distance between S 1 and S 2 : 

PK (S\S 2 )=31. 

Vector-like proximities allow more prominent description of dissimilarity between 
structures S 1 and S 2 : 

niS 1 ) = (MS 1 ),..., MS 1 ,..., MS 1 ) = (3,1,3,1,4,4,3,4,2), 



tt(S 2 ) = (MS 2 ), -,MS 2 , -,n 4 (S 2 ) = (2, 3, 2, 4, 3, 4, 1,3, 1), 
SJ(S\S 2 ) = (1,-2, 1,-3, 1,0, 2, 1,1), 
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V 

x(S\S 2 ) = 
y(S\S 2 ) = (y- 



Note, the coefficients — and 
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4 


o 


1 


-l 





-3 . -3 


-1 - 


-3 


-2 


-4 - 


-3 -3 


3 


-4 





1 


1 





-4 1 4 




-4 


-3 


-5 - 


-4 -4 


3 


4 




1 


-1 





-1 -2 -1 


3 - 


-1 




-1 


-1 -1 


1 4 -1 


5 


1 


1 




1 1 


3 


4 





1 


-1 





3 


4 





1 


-1 





(x- 3 ,x- 2 ,x- 1 


X , X , 


x 3 


) = (1,1,0,5,1,0), 


6 ,y-\y- A ,y- 3 




-l 


y\y 2 


,y 3 ,y i 


,y 5 ,y 6 ) 


(0,1,5,5,1,6,6,0,1,1,0,0). 







i(n-l) 



were not used in a; and y. 



5.4. Proximity for Trees. The following approaches to proximity (similarity/ 
dissimilarity) of trees are considered (e.g., [33 , [44 , [2 , 3 , [52 , [75 , [141 , [149 , 

[us], [nz]): 

(i) metrics (distance) (e.g., [33], [44], [149] . |157| ) including some special kinds 
of distances: alignment distance [50], isolated subtree distance [150] . top-down 
distance ( |141j . |167] h and bottom- up distance [156] ; 

(ii) tree edit distance (e.g., [31], [T4Q], [148], [154], [168], [169], [T70]); and 

(hi) common subtree, median tree or tree agreement, consensus (e.g., [2], [3], 
[52], [54], [75], [76], [142], [146], [162]). 

It is reasonable to note, trees are ordered structures. Thus, efficient (polynomial) 
computing algorithms have been suggested for building metrics/proximities of some 
kinds of trees (e.g., [49], [154], [155], [156]). 

Further, our simplified version of two-component proximity for rooted labelled 
trees is considered and used. Let T = {A',E') and T" = {A" ,E") be two rooted 
labelled trees (the root is the same in both trees) where A' and A" are the vertices, 
E' and E" are the arcs. Let us consider dominance parameter V (a, b) G (E* f]E"). 
The following three dominance cases can be examined: (i)) a — > b (a dominates b, 
i.e., a >- b), (ii) b — > a {b dominates a, i.e., a -< 6), and (hi) a and b are independent. 
Then dominance parameter is: 

{d 1 , if a -> 6, 
d 2 , if 6 a, 
d 3 , if a, 6 are independent. 

As a result, change parameter V(a, 6) € (i5' p| i?") can be define as the following: 

0, if rf( a ,6)e(i5' pi e") is not changed, 

1, if d(a,b)e(ET\E") is changed. 



P(a,b)e(E' f[E" 



Finally, the proximity of two trees is: 

p(T',T") = (p(A',A"),p(E / ,E")), 

where 

Pe (A',A") = l 



\A'f]A"\ 



p(E',E") = 



\A'\JA»\' 

J2(a,b)£(E> E") P(a,b)G(E> Q E") 

\(E>f)E")\ 



22 



MARK SH. LEVIN 



Evidently, the following properties are satisfied (Freshe axioms 1 and 2): 
(1)0 </>(£',£") <1, (2) p(E,E)=0. 

Note the described approach can be applied for digraphs as well. Now let us 
consider illustrative numerical examples for trees. 

Example 5.1: Trees T" and T" are presented in Fig. 5.2. Clearly, the general 
proximity is: p{T \T") = (0, 1). 



1 2 3 4 5 6 7 8 
T' • 



T" 



Fig. 5.2. Tree for example 1 



Example 5.2: Trees T' and T" are 
p(A',A") = 0. Table 5.2 and table 5.3 
for T 1 and T" (changes are depicted via 
p(T',T") = (0,JL). 



presented in Fig. 5.3. Here A' = A" and 
contain corresponding dominance factors 
'ovals'). Finally, the general proximity is: 
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Tabic 5.3. Dominance factor for T" 
1 2 3 4 5 6 7 




Example 5.3: Trees T and T" are presented in Fig. 5.4: A' = {1, 2, 3, 4, 5, 6, 7, 8}, 
A" = {1, 2, 3, 6, 8, 9, 10}, and A' f| A" = {1, 2, 3, 6, 8}. The proximity for vertex sets 
is: p{A',A") = \. 

Table 5.4 and table 5.5 contain corresponding dominance factors for A 1 1"| A" in 
T" and in T" (changes are depicted via 'ovals'). Finally, the general proximity is: 
p{T',T") = {\,\). 





Fig. 5.4. Tree for example 3 



Table 5.4. Dominance factor for T" 



Table 5.. Dominance factor for T" 
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1 






d 3 


© 


d 3 


3 








d 1 


d 1 


3 








d 1 


© 


6 








* 
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d 3 
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5.5. Proximity for Morphological Structures. Generally, morphological struc- 
tures are composite ones and it is reasonable to consider the corresponding their 
proximity as vector-like proximity. General structure (A) consists of the following 
parts: (i) tree-like system model T, (ii) set of leaf nodes P, (iii) sets of DAs for each 
leaf node D, (iv) DAs rankings (i.e., ordinal priorities) R, and (v) compatibility 
estimates between DAs I. Thus, the vector-like proximity for two morphological 
structures A Q , A^ can be examined as follows: 

~p (A", A") = (p*(T a ,T /3 ), p*(P«,P^), p*(D a ,D^), p'(R Q ,R0, P 1 ^^))- 
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Let us consider a simplified example: structure (A) is examined as a composition 
of two parts: (a) tree T = (A, E) (i.e., set of vertices and set of arcs), (b) rankings 
for each leaf node i ( U^iJj ). In our case, the proximity of two morphological 
structures is: 

p(A',A") = (p(A',A"),p(E',E"), Pr (A',A")). 

Here it is assumed sets of design alternatives for leaf vertices are not changed and 
tables of compatibility estimates are not considered. An example illustrate the 
approach. 

Example 5.4: Trees T and T" are presented in Fig. 5.5: A' = {1, 2, 3, 4, 5, 6, 7, 8}, 
A" = {1,2,3,4,5,6,7,9}, and A'\JA" = {1,2,3,4,5,6,7}. The proximity for ver- 
tex sets is: p(A',A") = |. Table 5.6 and Table 5.7 contain corresponding dominance 
factors for A' (J A" in T" and in T" (changes are depicted via 'ovals'). Finally, the 
proximity for trees is: p(T',T") = (§, ^-). 




A41 (2) Abi(1) Agi(2) A 7 i(3) Asi(3) 

A 42 (l) As 2 (3) A 62 (2) A 72 (2) A 82 (2) 

A 43 (3) A5 3 (2) As 3 (3) A 73 (l) As 3 (2) 

A 44 (2) A 54 (2) A 64 (l) A 74 (l) As 4 (l) 



Au(l) Asi(2) Asi(l) A n (3) A 9 i(l) 

A 42 (2) A 52 (3) Aa 2 (2) A 72 (2) A 92 (l) 

A 43 (3) Ab3(2) Aas (3) A 73 (2) A 93 (2) 

A 44 (3) A 54 (l) A 64 (2) A 74 (l) A 94 (2) 



Fig. 5.5. Trees with morphology for example 4 



Here the comparison process for rankings is based on rankings for common leaf 
vertices: L(T',T") = {4,5,6,7}. Thus, for each vertex above it is possible to 
use a metric/proximity. In our case (Fig. 14), the following metric (normalized 
Kendall- Tau distance) values are obtained: vertex 4: <5 4 (A',A") = -|, vertex 5: 



S 6 (A',A") 



vertex 6: <5 6 (A',A") 



and vertex 7: <5 7 (A',A") 



The 



resultant general proximity for rankings can be computer as an average value: 



MA', A") 



ieL(T',T" 



Si(A',A") 



\L(T',T")\ 



0.33. 



Finally, the general proximity is: 



p(A',A") = (p(A',A"),p(E',E"),p r (A,A")) = (0.22,0.1,0.33). 



Note, proximity of compatibility estimates can be added into the proximity vector 
above as well. 
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Table 5.6. Dominance factor for T' 
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Table 5.7. Dominance factor for T" 
1 2 3 4 5 6 7 
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6. Underlying Problems 

6.1. Multicriteria Ranking. Let A — {Ai, Aj, A n } be a set of alterna- 
tives and C = {Ci, C p , C m } be a set of criteria. An estimate vector = 
{zji, Zj p , Zjm} a s a result of an assessment upon criteria above is given for 
each alternative Aj. In addition, criteria weights {^ p } are used as well. The prob- 
lem consists in comparison of the alternatives and forming for each alternative an 
ordinal quality estimate (priority). This problem belongs to a class of ill-structured 
problems by classification of H. Simon |144) . In the paper, a modification of out- 
ranking technique Electre [139] has been used (the modification was suggested in 

ESI). 

6.2. Knapsack Problem. The basic problem is (e.g., [57], [86j): 

m 

max^ CiXi 

i=\ 

m 

s.t. aiXi < b, Xi G {0, 1}, i = 1, m 

i=l 

and additional resource constraints J^iLi a i,kXi < k = 1,1; where X{ = 1 if item 
z is selected, Cj is a value ("utility") for item i, and is a weight (or required 
resource). Often nonnegative coefficients are assumed. The problem is NP-hard 
[57) and can be solved by the following methods: (i) enumerative methods (e.g., 
Branch-and-Bound, dynamic programming), (ii) approximate schemes with a lim- 
ited relative error (e.g., [86]), and (iii) heuristics. 
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A basic version of knapsack problem with minimization of the objective function 

is: 

m 

min y ] CiXi 
i=i 



s.t. a i x i >b, Xi 6 {0, 1}, i = 1, m. 

i=l 

For multiple criteria statements it is reasonable to search for Pareto-efficient 
solutions. Mainly, heuristics, evolutionary algorithms, dynamic programming, and 
local search methods are applied to multicriteria knapsack problems (e.g., [15] . 
[48] . [SB] . [92] . [124] . |172| ). A recent survey on multicriteria knapsack problem is 
contained in [124] . 

6.3. Knapsack Problem and Compatibility. Now let us consider item depen- 
dence in knapsack problem [101] . Here the following is considered: binary compati- 
bility of items as a symmetric binary relation (i.e., the selected subset has to contain 
only compatible items, 1 corresponds to compatibility of items and correspond to 
incompatible case). Thus, the consideration of compatibility for knapsack problem 
leads to searching for a clique ("profit clique") with the following properties: 

(i) maximum total profit of the selected items: niax^™ =1 a, 

(ii) restricted total weight of the selected items: Y^i=i a i — 

Table 6.1 and Table 6.2 contain descriptions of numerical examples for the fol- 
lowing problems (this example is a modification of the example from [101] ): 



Table 6.1. Initial data and results 



i 


Ci 




Knapsack 


"Proht 


Maximal 








problem 


clique" 


clique 


1 


5 


2 








2 


12 


5 








3 


4 


2 








4 


2 


1 








5 


3 


2 


•k 






6 


4 


3 








7 


1 


2 








8 


1 


3 








C 




26 


23 




A 




12 


11 




M 




5 


4 


6 



(a) basic knapsack problem, 

(b) "profit clique" (i.e., knapsack problem and compatibility of the selected 
items), and 

(c) maximal clique (maximizing the number of selected compatible items). 
The following notations are used: 

(1) M is the number of items in a solution (i.e., the number of the selected 
items); 

(2) b = 13 is the right-side constraint (additional constraints are not applied); 

(3) symbol * corresponds to a selected item; and 

(4) characteristics the solution are: C — YhLi c i x ii A — Eti a i x i- 
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Clearly, the considered problems are NP-hard. Thus, heuristics are used. Mul- 
ticriteria versions of "profit clique" problem can be examined and used as well. 



Table 6.2. Compatibility 



i 


1 


2 


3 


4 


5 


6 


7 


8 


1 




1 





1 





1 


1 


1 


2 


i 







1 


1 


1 


1 


1 


3 


l 


b 




1 


1 


1 


1 


1 


4 


l 


i 


i 




1 


1 


1 


1 


5 








l 







1 


1 





6 


l 


i 





1 







1 


1 


7 


l 


i 


l 


1 


1 


i 




1 


8 


l 


i 





1 





l 


i 





6.4. Multiple Choice Knapsack Problem. In the case of a multiple choice prob- 
lem, the items (e.g., actions) are divided into groups and we select elements from 
each group while taking into account a total resource constraint (or constraints) 
(e.g., [57], [73], [86], [132]): 

m qi 

CijXij 

i=l j=l 

m qi qi 

s.t. CijXij <b,*^2 x ij = 1) i = l,m, Xij € {0, 1}. 

i=l 3 = 1 j=l 

In the case of multicriteria description, each element (i.e., has a vector 

profit : Cjj = (cjj , .. . , cf j , ... , c£ j ) . A version of multicriteria multiple choice problem 
was presented in ([110]', QT8], [120] ): 

m qi 

max X! X! '{<■'' :i > Vf = 1, r 

«=1 3=1 

si. / ^ / J ajj^jj < 6, J^^ii = 1, i = l,m, G {0, 1}. 

1=1 3=1 3=1 

Evidently, in this case it is reasonable to search for Pareto-efficient solutions (by 
the vector objective function above). Here the following solving schemes can be 
used f |118j . [120] ~1 : (i) dynamic programming, (ii) heuristic based on preliminary 
multicriteria ranking of elements to get their priorities and step-by-step packing 
the knapsack (i.e., greedy approach), (iii) multicriteria ranking of elements to get 
their ordinal priorities and usage of approximate solving scheme (as for knapsack 
problem) based on discrete space of system excellence (i.e., quality lattice as in 
HMMD [101] ). In the article, greedy heuristic above is used later. 

6.5. Morphological Design. A brief description of HMMD is a typical one as 
follows (e.g., |101| . [105] . [106] . |110] V The composite (modular, decomposable) 
system under examination consists of the components and their interconnections 
or compatibilities. Basic assumptions of HMMD are the following: (a) a tree-like 
structure of the system; (b) a composite estimate for system quality that integrates 
components (subsystems, parts) qualities and qualities of interconnections (here- 
inafter referred as 'IC') across subsystems; (c) monotonic criteria for the system 
and its components; and (d) quality of system components and IC are evaluated on 
the basis of coordinated ordinal scales. The designations are: (1) design alternatives 
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(DAs) for nodes of the model; (2) priorities of DAs (r = 1, k; 1 corresponds to the 
best level); (3) ordinal compatibility estimates for each pair of DAs (w = 0,1; I 
corresponds to the best level). The basic phases of HMMD are: 1. design of the 
tree- like system model (a preliminary phase); 2. generating DAs for model's leaf 
nodes; 3. hierarchical selection and composing of DAs into composite DAs for the 
corresponding higher level of the system hierarchy (morphological clique problem); 
and 4- analysis and improvement of the resultant composite DAs (decisions). Let S 
be a system consisting of m parts (components): -P(l), P(i), • P(m). A set of 
design alternatives is generated for each system part above. The problem is: 

Find a composite design alternative S = S(l) * ... * S(i) * ... * S(m) of DAs 
(one representative design alternative S(i) for each system component/part P(i), 
i = l,m) with non-zero IC estimates between design alternatives. 

A discrete space of the system excellence on the basis of the following vector is 
used: N(S) = (w(S);n(S)), where w(S) is the minimum of pairwise compatibility 
between DAs which correspond to different system components (i.e., V Pj 1 and 
Pj 2 , 1 < ji 7^ 32 < Tn) in S, n(S) = (ni, n r , ...rife), where n r is the number of 
DAs of the rth quality in S (X^ r=1 n r = m). As a result, we search for composite 
system decisions which are nondominated by N(S). Here an enumerative solving 
scheme (e.g., dynamic programming) is used (usually m < 6). 

Fig. 6.1 and Fig. 6.2 illustrate the composition problem (by a numerical example 
for a system consisting of three parts S = X * Y * Z). Priorities of DAs are shown 
in Fig. 6.1 in parentheses and are depicted in Fig. 6.2; compatibility estimates are 
pointed out in Fig. 6.2). In the example, the resultant composite decision is (Fig. 
6.1, Fig. 6.2): Si = A 4 * Y 2 * Z 2 , JV^) = (2; 1, 1, 1). 



1 



S =X*Y*Z 

Si = X\ -k Yi -k Z 2 



X Y | 



O X 3 (1) ®Y 3 {3) 1<S>Z 3 (3) 

k) (xM} k)Y 4 (l) 

Fig. 6.1. Example of composition 



Z 2 



Xi) 



x 1 



x 2 
x 3 



Y 2 
5ht 



Y 1 Y 3 



Fig. 6.2. Concentric presentation 



7. Median/Consensus Problems and Aggregation Problems 

7.1. Sets. Evidently, sets are the basic structures. Let us consider the case of m 
sets. The extended median/consensus for m sets {A\, ...,Ai,..., A m } is the follow- 
ing. Analogically, Ve G 1J™ 1 Aj) two attributes are examined: "profit" / "utility" c e > 

0, required resource b e > 0. Let R[Ai} = ${Ai} Q fXiLi (or R{Ai} = ^{Ai}) De 
a basic "consensus" set and ^2 e€RAB < b, where b is a total resource. The problem 
of "extended median/consensus" is: 



max °e 
eewc((U™ =1 AA\ R{Ai) ) 
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S.t. be < b. 

e£WC((\J™ =1 Ai)\R {A . } ) 

This model corresponds to basic knapsack problem. 

Further, aggregation problems for sets will be examined. First, let us consider 
two-set case for sets A, B. Here the following proximity/metric is used: p(A, B) > 
0, p(A, B) = p(B, A). For example, the following simple metric can be used: p(A, B) = 
\Af]B\/\A[jB\. 

A median-like subset (median-like consensus model by Kendall etc. [90]) is: 
M AB = arg min (p(M, A) + p(M, B)). 

{M} 

The case of the extended median/consensus for two sets A , B is similar. Here 
two attributes are examined Ve G (A{JB): "profit" / "utility" c e > 0, required 
resource b e > 0. Let Rab — Sab (or R AB — Mab) be a basic "consensus" set 
and ^2e£R A B — ^' wnere b is a total resource. The problem of building the "extended 
median/consensus" is: 

max c e 

eeWCaA(JB)\R AB ) 

s.t. Y be < b. 

eeWC({A\JB)\R AB ) 

This model corresponds to basic knapsack problem. 

Now let us consider the multi-set case. The extended median/consensus for m 
sets {Ai, Ai, A m } is the following. Analogically, Ve 6 IJfci two attributes 
are examined: "profit" / "utility" c e > 0, required resource b e > 0. Let R{A t } = 

S{Ai} Q Hi™ i Ai (or R{Ai} = M{Ai}) be a basic "consensus" set and J2 e eR AB ^ e — 
b, where b is a total resource. The problem of "extended median/consensus" is: 

max c e 

eGWC(({JT =1 Ai)\R {At} ) 

s.t. Y b e < b. 

eWC((Ur =1 Ai)\R {A . } ) 

This model corresponds to basic knapsack problem. 

7.2. Rankings. Median/consensus of rankings is one of the basic problems in de- 
cision making (e.g., QjZ], [3S], [55], [55], [H], [55J, [50]). Here three methods are 
briefly described: (i) median consensus method based on assignment problem (e.g., 
[34], [35], [36], [95]); (ii) heuristic approach (e.g., [16], [95], [152]); and (iii) method 
based on multiple choice problem (e.g., |100| . [101] ). 

Method 1. The median consensus method based on distance (usually: Kendall- 
Tau distance) and assignment problem has been studied by Cook et al. (e.g., 
[34], [35], [36]). Our version of the approach (for layered sets) is the following. 
Let A — {1, i, n} be the initial set of elements (alternatives, objects). The 
number of layers equals m {k = 1, m). There are p initial rankings of set A: S 1 ,..., 
S x ,..., S>*. Thus, S x = UfeLi A k- Let (i G A) be the priority of i in S x , i.e., the 
number of corresponding layer: r x — k if i G A^. 
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The resultant ranking (consensus) is S a = {A\, A%, A^}, and correspond- 
ing consensus priorities are rf Vi £ A. The following binary variable will be used: 

f 1, if rf = k or i e A%, 
tk y 0, otherwise. 

The assignment problem for finding the consensus is (case of layered set): 

n rn fi 

i=l k=l \=1 

m 

s.t. ^2 x ik = 1, i = 1, n; Xik £ {0, 1}. 

k=l 

Generally, polynomial algorithms exist for basic assignment problem (e.g., [57], 
[96] V The obtained version of assignment problem is more simple and can be 
solved by an evident greedy algorithm: selection of the closest layer Mi e A. Let us 
consider a numerical example (Table 7.1): 



Table 7.1. Illustrative example for rankings 



i e A 


r\ (S 1 ) 


rf (S 2 ) 


rf (S 1 ) 


rf (S a ) 




1 


3 


3 


3 


3 


3 


2 


1 


1 


1 


1 


1 


3 


3 


1 


2 


2 


2 


4 


1 


2 


1 


1 


1 


5 


4 


4 


3 


4 


4 


6 


4 


4 


4 


4 


4 


7 


3 


3 


4 


3 


3 


8 


4 


4 


4 


4 


4 


9 


2 


2 


2 


2 


2 



(a) initial set of elements A — {1, 2, 3, 4, 5, 6, 7, 8, 9}, 

(b) three rankings (four layer are examined): 

5 1 = {A\ = {2, 4}, A\ = {9}, A\ = {1, 3, 7}, A\ = {5, 6, 8}}, 

5 2 = {A\ = {2, 3}, A\ = {4, 9}, A\ = {1, 7}, A\ = {5, 6, 8}}, 

5 3 = {A\ = {2, 4}, A\ = {3, 9}, ^ = {1, 5}, A\ = {6, 7, 8}}. 

The resultant ranking based on assignment problem above is (Table 10): 
S a = {A* = {2, 4}, A% = {3, 9}, A% = {1, 7}, A£ = {5, 6, 8}}. 

Method 2. A basic heuristic approach has been suggested in [16]. The method 
is widely used (e.g., |152j ). Let us consider a simplified version of heuristic to find 
the corresponding solution S a : 

= / KELi r A )M if (ELi r, A )A* - [(Ea=i < 0.5, 

1 I KESU^VmI+I. otherwise. 

Thus, solution is (Table 10): 

S» = {At = {2, 4}, A% = {3, 9}, A% = {1, 7}, A£ - {5, 6, 8}}. 
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Method 3. Now a more complicated aggregation process based on multiple choice 
problem is examined (initial rankings {>S' A |A = 1,/i} are mapped into a fuzzy ag- 
gregated ranking S a ) {S x } S a [TUT] : 

h{S a | s a G {S | 77(5, 5 A ) r< T] VA = 1, fi} — > max, 

where is an attribute (quality) of the resultant ('overage') structure S a ; r\ Q is a 
vector- like proximity. The problem is depicted in Fig. 7.1. 




Fig. 7.1. Aggregation [101J 

The following notations are used: an is the number of initial structures, in which 
element !g4i,I=l,m (layers); vector = (^i, £j m ) defines frequencies 

of belonging of element i to layers {Ai, ...,Ai, ...,A m }, where £y = 2ii (it is a 
membership function of element i to layer I = 1, .., m). Let us denote 5^ as a set 
of intervals {di}. The resultant quality of SJ is based on the following entropy-like 
function: 

n n 1 

1 = 1 4=1 ' 1 

Next, a modular vector as the proximity is used: z = (z 1 , z k , 0, 0). Finally, 
the problem is: 

n 

^i^(S Q ) — ■» max, z(S A ,S Q ) r< z Q ,VA = 1, 

i=l 

With respect to zero- valued components z fc+1 , ...,z fe+ , it is possible to define a set 
of admissible variants to intervals {dig \ 9 — 1, ■ ••,%)• Thus, we reduce the model 
to the following modification of multiple-choice problem ([ST], [86]): 

n qi 

^^Hie{die)K ig — ► max, 

i=l 0=1 

n qi qi 

EX/X, < 22^, P = h—,k, }] K t e = 1, i = l,...,n,n ie G {0,1}, 

r>p i=l 9=1 r =P 0=1 

where b r iS is the sum of components £j, which are differed from d\ g {df g ) by r. A 
version of the described aggregation scheme has been implemented in DSS COMBI 

(HE], EH). 

7.3. Trees. Here initial information consists in a set of trees. Usually four basic 
approaches are considered: 

(1) maximum common subtree (e.g., [2], [54]); 

(2) median/agreement tree (e.g., [3], [JJ], E>2], E3, [IS], [133], Q35J); 

(3) compatible tree (e.g., [T5], [54], [57]): and 

(4) maximum agreement forest (e.g., [30], [55], [138] , [164] , |165] V 
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Mainly, the problems above correspond to class of NP-hard problems (e.g., [67], 
[133] ). As a result, heuristics, approximation schemes, and enumerative methods 
are used. 

Thus, the aggregation problem for set of initial trees {T} = {T 1 , ...,T l , ...,T m } 
can be considered as follows (addition strategy I): 

Stage 1. Searching for a median-like tree (i.e., "kernel"): 

m 

T a °° = arg mn£p(r,f)). 

* ' i=i 

Stage 2. Generation of a set of additional elements (nodes and/or edges). 
Stage 3. Addition of elements to T a " (knapsack- like problem). 

Evidently, in the case of vector- like proximity p(T', T"), T a " has to be searched 
for as Pareto-efficient solution(s). On the other hand, it is reasonable to consider 
some heuristic algorithms for building the "kernel" , for example: 

m— 1 

K = U {T l [\T l=1 )- 

i=l 

7.4. Morphological Structures. The basic aggregation problem for set of initial 
structures {A} = {Al, A 1 , A m } can be considered as follows (addition strategy 
I): 

Stage 1. Searching for a median-like tree (i.e., "kernel"): 

m 

A Q ^ = arg mm£>(A,A*)). 

i— 1 

Stage 2. Generation of a set of additional elements (nodes and/or edges). 
Stage 3. Addition of elements to A a " (knapsack- like problem) . 

Evidently, in the case of vector-like proximity p(A', A"), A a " has to be searched 
for as Pareto-efficient solution(s). Generally, morphological structures (morpholog- 
ical structures with compatibility) are very complicated composite structures: A= ( 
T,P,D,R,I ). Here, it may be reasonable to consider the following heuristic solving 
scheme: 

Stage 1. Aggregation of sets of systems parts {P}. 
Stage 2. Aggregation of sets of DAs {D}. 
Stage 3. Aggregation of sets of compatibility estimates {I}. 
Stage 4- Aggregation of tree-like models {T}. 



8. Illustrative Applied Numerical Examples 
The list of examined applied examples is presented in Table 8.1. 
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Table 8.1. List of examined applied examples 



Applied example 


Aggregation strategy 


Taking into 


References 






account 


for 






compatibility 


example 


1. On-board tele- 


addition (strategy 1) 


None 


11 1 rtl 


metry system 








2. Notebook 


combined strategy 


None 












3. Zig-Bee protocol 


(i) addition (strategy 1) 


None 


liijj 




(ii) deletion (strategy II) 


None 




4. Integrated se- 


addition (strategy I) 


None 


[117] 


curity system 








5. Plan of students 


"kernel" -based strategy 


None 


h n1 l 


art activity 








6. Combinatorial 


"kernel" -based strategy 


None 


Lioi] 


investment 








7. Common educa- 


(i) addition (strategy I) 


None 




tional course 


(ii) median-based strategy 


None 




8. Educational 


extended strategy IV 


Yes 


m, mi 


course on design 






[T02]. [106] 


9. Configuration of car 


median-like strategy 




una 


10. Applied Web-based 


(i) kernel-based strategy 


Yes 


UJL4J 


information system 


(ii) design strategy (IV) 







8.1. ZigBee Communication Protocol. Here an applied example for aggrega- 
tion of ZigBee protocol solutions from |119j is briefly described. A general frame- 
work for aggregation of two solutions is depicted in Fig. 8.1. A brief description of 
Zigbee protocol components is the following [119] : 

1. Interference avoidance A: 1.1. Startup Procedure of Channel Acquisition X . 

1.2. Channel Hopping J. 

2. Automated/distributed address management B. 

3. Group addressing /. 

4- Centralized data collection C: 4.1. Low-overhead data collection by ZigBee 
Coordinator G. 4.2 Low-overhead data collection by other devices H. 4-3- Many- 
to-one routing Q. 4-4- 6L0WPAN multicast/broadcast support V. 4-5- Source 
routing P. 

5. Network scalability D. 

6. Message size E. 

7. Standardized commissioning Z . 

8. Robust mesh networking: 6L0WPAN approach U . 

9. Cluster Library support L. 

10. Web services support W. 

The above-mentioned system components have implementation alternatives (e.g., 
B 1 ,B 2 ). 

The following initial solutions are considered: 
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(a) an expert-based forecast 1 (S4 from |119j ) (Fig. 8.2), 

(b) forecast based on using knapsack problem 2 (<£> from |119j ) (Fig 



8.3). 



Expert-based 
solution 1 



Substructure 









Aggregation 
strategy I 
(extension) 



e 1 




Knapsack-based 
solution 2 




Aggregation 
strategy II 
(deletion) 




Fig 8.1. Example: aggregation for ZigBee protocol 



Expert-based forecast 6 1 : ZigBee/IP (6L0WPAN) 2010 S A 

La IbIiL 



B 1 h 



X U 

X! Jl 



c 



Id Ie lz [u II \w 

D 2 E 3 Z x U x Lx Wx 

u 2 



G \h\Q \v 

Gx Hx Qi Vx 

V 2 

Fig. 8.2. Structure of direct expert-based forecast O 1 |119) 



Knapsack-based forecast 2 : $ 



7a Ib ii lc 

B 2 h 

X 
Xx 



Id Ie lz 

D2 E2 Z\ 

Qi Pi 

Fig. 8.3. Forecast based on knapsack problem 2 [119] 



The corresponding substructure of protocol and superstructure of protocol 
are shown in Fig. 8.4 and in Fig. 8.5. 
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Substructure 
B 2 h 



A 



C 



X 
Xi 



~\q 

Qi 



Id [e \z 

D 2 E 2 Zi 



Fig. 8.4. Substructure of forecasts 6 



Superstructure 



.4 



X I J B2 

X! Jl 



c 



Id lz lu II \w 



g\h\q\ 

Gi m 



D 2 

~v~[p 



E 2 
E, 



Zi 



Ui 

U 2 



L x Wi 



Y 2 



Pi 
P 2 



Fig. 8.5. Superstructure of forecasts 

Further, let us consider aggregation processes. A list of addition operations 
(for strategy I) is presented in Table 8.2 (designation of operation from [119] is 
depicted in parentheses). Here and hereafter the following attributes (criteria) 
for an assessment of the addition operations (deletion operations) are used |119j : 
(1) cost Ti; (2) required time for implementation T2; (3) performance T3; (4) 
decreasing a cost of maintenance T 4 ; (5) scalability T 5 ; (6) reliability T 6 ; (7) 
mobility T7; and (8) usability value Tg. An ordinal scale [1,5] is used for each 
criterion: 1 corresponds to "strong negative effect", 2 corresponds to "negative 
effect", 3 corresponds to "no changes", 4 corresponds to "positive effect", and 5 
corresponds to "strong positive effect" . Priorities {r^} are obtained via multicriteria 
ranking (Electre-like method |119j ). 

Table 8.2. Addition operations (estimates, priorities) 



i 


Improvement 


Vari- 


Ti 


T 2 


T 3 


T 4 


T 5 


T 6 


Ty 


T 8 


Priori- 




operation 


able 


















ties r. 


1 


Ji ($10) 


Xi 


2 


3 


4 


4 


3 


4 


4 


3 


1 


2 


Bi ($13) 


x 2 


3 


3 


3 


4 


3 


3 


4 


3 


3 


3 


J7i&J7 2 ($17) 


%3 


3 


3 


4 


3 


4 


3 


3 


3 


3 


4 


Li ($3) 


X4 


3 


4 


3 


4 


3 


5 


3 


4 


1 


5 


Wi ($15) 


X5 


3 


3 


2 


3 


4 


3 


3 


5 


2 



Thus, the addition problem (simplified knapsack problem) is: 



max > CiXi 

i=l 



5 

s.t. diXi < b, Xi G {0, 1}. 

i=l 
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Cost estimates are (by criterion Ti) used as {a^}, priorities {r^} are used as (trans- 
form to) {ci}, and b = 8.00. 

A resultant solution for strategy I is depicted in Fig. 8.6 (xx = 1, X2 = 1, 
X3 = 0, X4 = 0, £5 = 1). Here compatibility estimates between design alternatives 
for system components are not considered. 



J (strategy I) 



.4 



Bi h 
B 2 



C 



Id Ie lz \w 

Zx Wi 



Do E 2 



[x u " 2 W]p 

Xx Ji Qi Px 

Fig. 8.6. Aggregated solution 9 7 based on strategy I 

A list of deletion operations (for strategy II) is presented in Table 8.3 (designation 
of operation from |119j is depicted in parentheses). 

Table 8.3. Deletion operations (estimates, priorities) 



i 


Deletion 


Vari- 


Ti 


T 2 


T 3 


T 4 


T 5 


T 6 


T 7 


T 8 


Priori- 




operation 


able 


















ties Ti 


1 


Ji ($10) 


Xi 


2 


3 


4 


4 


3 


4 


4 


3 


1 


2 


Bx ($13) 


X2 


3 


3 


3 


4 


3 


3 


4 


3 


3 


3 


Qx (*7) 


X3 


3 


4 


3 


3 


3 


2 


3 


3 


4 


4 


Pi (* 8 ) 


X4 


3 


4 


4 


4 


3 


4 


3 


3 


2 


5 


E 3 ($ 14 ) 


X5 


3 


3 


3 


3 


5 


3 


3 


3 


2 


6 


Li ($3) 


X 6 


3 


4 


3 


4 


3 


5 


3 


4 


1 


7 


U1&U2 ($i 7 ) 


X 7 


3 


3 


4 


3 


4 


3 


3 


3 


3 



Thus, the deletion problem (knapsack problem with minimization of the objec- 
tive function) is: 

6 

min^^ axi 

i=X 



s.t. ciiXi > b, Xi £ {0, 1}. 



i=X 



Cost estimates are (by criterion Ti) used as {a^}, priorities {r^} are used as (trans- 
form to) {ci}, and b = 8.00. 



G n (strategy II) 



B 2 h 



Id Ie lz 

D 2 E 2 Z\ 

V E 3 



IlIw 



A 



X ] J 

Xx Jx 



C 



G [h 1 

G x Hx 1 



Fig. 8.7. Aggregated solution Q 11 based on strategy II 
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A resultant solution based on strategy II is depicted in Fig. 8.7 [x\ = 0, x% = 1, 
2:3 = 1, Xi = 1, X5 = 0, xe = 0, X7 = 1). Here compatibility estimates between 
design alternatives for system components are not considered. 

8.2. On-Board Telemetry System. Here a numerical example for on-board 
telemetry system is considered from |116| . The initial morphological structure for 
on-board equipment is the following (Fig. 8.8): 
1. On-board equipment S = D * E * P. 

1.1. Power supply D — X *Y * Z: 1.1.1. stabilizer X: X x (standard), X 2 
(transistorized), A 3 (high-stability); 1.1.2. main source Y: Y\ (Li-ion), Y% (Cd- 
Mn), Y 3 (Li); 1.1.3. emergency cell Z: Z x (Li-ion), Z 2 (Cd-Mn), Z 3 (Li). 

1.2. Sensor element E = O * G: 1.2.1. acceleration sensors /: I\ (ADXL), I 2 
(ADIS), I 3 (MMA); 1.2.2. position sensors O: O x (SS12), 2 (SS16), 3 (SS19), 
O4 (SS49), O5 (SS59), 6 (SS94); 1.2.3. global positioning system (GPS) G: G 1 
(EB), G 2 (GT), G 3 (LS), G 4 (ZX). 

1.3. Data processing system P = H * C * W: 1.3.1. data storage unit H: Hi 
(SRAM), H 2 (DRAM), H 3 (FRAM); 1.3.2. processing unit (CPU) C: d (AVR), 
C 2 (ARM), G 3 (ADSP), G 4 (BM); 1.3.3. data write unit W: W x (built-in ADC), 
W 2 (external ADC I2C), W 3 (external ADC SPI), W 4 (external ADC 2W), W 5 
(external ADC UART(l)). 



D = X*Y*Z 



) D l =X 2 *Y 2 *Z 2 
D 2 =X 3 *Y 3 * Z 3 



O 



X Y 



Z 



> S(A) = D*E*P 
E = I-kO-kG 



P = H-kC-kW 



E 1 =h*0 1 -k G 4 
E 2 = Ii*O s * G 4 
E 3 = h * Ox * G 4 



X x {2) Y x {2) Z x {2) 

a 2 (i) y 2 (i) z 2 (i) 

A 3 (3) Y 3 {2) Z 3 (3) 



6 
h{2) 

Hi) 





= H 2 *C X *W 2 


Pi 


= H 2 * Gi * W 5 


P3 


= H 3 * Gi * W 2 


Pi 


= H 3 *C X *W 5 


i H 





°6 G < 

Oi(l) G x (l) H x (3) Gr(l) Wi(3) 

2 (2) G 2 (2) tf 2 (l) G 2 (2) W 2 (l) 

3 (2) G 3 (3) iT 3 (2) G 3 (3) W 3 (2) 

4 (3) G 4 (l) G 4 (3) W 4 (3) 

O s (l) W 5 (l) 
Oe(3) 

Fig. 8.8. Basic hierarchy for on-board telemetry system |116j 

In |116j . 24 resultant solutions have been obtained via HMMD (A x , A 24 ). In 
the example, 6 initial solutions are examined (in parentheses designation from [116] 
in pointed out): 

S X {A X ) = D X *E X *P X = (X 2 * Y 2 * Z 2 ) * (/1 * Qi * G 4 ) * (fT a *C X * W 2 ), 
S 2 (A 24 ) =D 2 *E 3 *P 4 = (X 3 * Y 3 * Z 3 ) * (I 3 * Q x * G 4 ) * (ff 3 * G 4 * W 5 ), 
S 3 (A n ) =D 1 *E 3 *P 3 = (X 2 * Y 2 * Z 2 ) * (J 3 * Q x * G 4 ) * (ff 3 * G 4 * W 2 ), 
S 4 (A 19 ) =D 2 *E 2 *P 3 = (X 3 *Y 3 * Z 3 ) * (Jr * Q 5 * G 4 ) * (ff 3 * C x * W 2 ), 
^5(^23) =D 2 *E 3 *P 3 = (X 3 * Y 3 * Z 3 ) * (J 3 * Qx * G 4 ) * (H 3 * G x * W 2 ), 
S 6 (A 12 ) =D X *E 3 *P 4 = (X 2 * Y 2 * Z 2 ) * (J 3 * Qx * G 4 ) * (Jf 3 * d * W 5 ). 
In Fig. 8.9 and Fig. 8.10, supersolution and subsolution are depicted. 
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X Y 

X2MY2 ^ 

X3 Y 3 



Z 

z 3 



I Q G H C W 



ft 



ft 



1 Hi* 
H 3 



ft 



f w 2 
w 5 



Fig. 8.9. Supersolution 



X Y 



Z 



I Q G H C W 



G4 Ol 



Fig. 8.10. Subsolution 8 

The obtained subsolution contains only two elements. Thus, the design of "sys- 
tem kernel" is based on the special method with a = 0.6 (Fig. 8.11). 



X Y 



z 



I Q G H C W 



ft 



ft 



ft 



Fig. 8.11. "System kernel" K 
Further, the addition strategy as simplified multiple choice problem is used: 



3 2 

CijXij 

i=i 3=1 



3 2 2 

s.t. 2J <HjXij < b,}^ Xij = lVt = 1, 3; Xij e {0, 1}. 

i=l j=l j=l 

A list of addition operations is presented in Table 8.4 |116j . Here cost estimates are 
based on expert judgment, priorities {r^} were computed via Electre-like method 
in [116] and transformed to {cj}, b = 9.00. 

A resultant solution is depicted in Fig. 8.12 (in = 1, X12 — 0, X21 — 1, X22 = 0, 
.T31 = 1, S32 = 0). Note, compatibility estimates between design alternatives for 
system components X,Y,Z are not considered. 



Table 8.4. Addition operations 



Improvement 


Binary 


Cost 


Priorities 


operation 


variable 




n 


x 2 


x n 


3 


1 


A 3 


xu 


4 


3 


Y 2 


X21 


2 


1 


Y 3 


X22 


3 


2 


Z 2 


X31 


2 


1 


z 3 


X32 


3 


2 
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X Y Z I 



G H C W 



x 2 ^ 



Ci w 2 



Fig. 8.12. Resultant aggregated solution 



8.3. Continuation of Example for Notebook. Let us examine the final part of 
the example for notebook. Here the combined aggregation strategy (strategy III) 
is considered for two cases: 

(1) "system kernel" as an extension of substructure: 

K' = B 1 * Ri * V 3 * Ji ★ Ei * Oi * Di * Ai * G\ * Lx * Q 2 ; 

(2) "system kernel" K* based on multicriteria selection of the "best" design 
alternatives for each system component. 

Let us consider case 1. Here the aggregation strategy as modification of "system 
kernel" K' can be applied. A set of candidate modification operations are the 
following: 

1. addition operations: 1.1. addition for U: U\ or U 2 or f7 3 , 1.2. addition for 
F: F x or F 2 , 1.3. addition for P: P 2 or P 3 or P 4 ; 

2. correction operations: 2.1. replacement Bx => B 2 , 2.2. replacement V3 =>■ V4, 
2.3. replacement Ax =>• A3. 

Table 8.5 contains the list of modification operations above, their estimates (or- 
dinal expert judgment) and corresponding binary variables. 



Table 8.5. Modification operations 



Operations 


Binary 


Cost 


Priorities 




variable 


a ij 


r ij 


1. Addition 








1.1. Ux 


Xxx 


3 


2 


1.1. u 2 


X12 


2 


3 


1.1. U 3 


X13 


4 


1 


1.2. Fx 


X 2 X 


2 


2 


1.2. F 2 


x 22 


3 


1 


1.3. P 2 


X31 


3 


1 


1.3. P 3 


X32 


2 


1 


1.3. P 4 


X33 





2 


2. Replacement 








2.1. Pi => P 3 


X41 


4 


1 


2.1. None 


X42 





2 


2.2. v 3 => y 4 


X51 


3 


1 


2.2. None 


X52 





2 


2.3. Ai A 3 


xei 


2 


1 


2.3. None 


Xf>2 





2 



The following simplified multiple choice problem is used (dj = 3 — r^, b — 11.00): 

6 9i 
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6 <?» g, 

s.t. aijXij < b, = IVi = 1, 6; Xij G {0, 1}. 

i=i j=i j=i 

A resultant computer solution S is depicted in Fig. 8.13 (in = 1, x\2 = 0, 
xi 3 = 0, X21 = 0, a; 2 2 = 1, x 3 i = 0, x 32 = 1, £33 = 0, Z41 = 0, a; 42 = 1, ar 5 i = 0, 
X52 = 1, xqi = 1, xq2 = 0). Here a greedy algorithm was used. Note, compatibility 
estimates between design alternatives for system components are not considered. 




Fig. 8.13. Solution S lc based on modification of "system kernel" K' 



Now let us consider case 2. Here building of "system kernel" is based on multi- 
criteria selection and/or expert judgment. The basic structure of "system kernel" 
is: B, U, R, V, O, F, D, G. For each system component above, it is possible to 
consider a selection procedure to choose the "best" system element (while taking 
into account elements of the initial solution or additional elements as well). 

Table 8.6 contains design alternatives for the selected components of "system 
kernel" structure above including ordinal estimates (expert judgment, the small- 
est estimates correspond to the best situation) and the resultant priorities. The 
following criteria were used: cost (Ti), usefulness (T2), experience (T3), prospec- 
tive features (T4). As a result, the following "system kernel" K* is obtained: 
K* = B 2 * U 2 * R 3 * V 3 * E 1 * Ox * F 2 * D x * G x . 

Further, the system correction process is based on the following operations: 

1. addition: 1.1. Ax, 1.2. Pi, 1.3. Lx\ 

2. deletion: 2.1. Ex] 

3. replacement: 3.1. B 2 B 3 , 3.2. U 2 => U 1 , 3.3. Ox 3 . 

Table 8.7 contains the list of modification operations above, their estimates (or- 
dinal expert judgment) and corresponding binary variables. 

The following multiple choice problem is used (cy = 3 — ry, b = 9.00): 

7 2 

i=X j = X 
7 2 2 

s.t. a^Xij <b,*^2 Xi i = = 7 ! x v e {°> !}• 

i=i j=i j=i 

A resultant computed solution S 2c is depicted in Fig. 8.14 (xxx = 1, 2:12 = 0, 
x 2 i = 0, .t 22 = 1, x 31 = 1, x 32 = 0, xax = 1, 2:42 = 0, ar 5 i = 1, x 52 = 0, a; 6 i = 0, 
a^62 = 1, X71 = 1, £72 = 0). Here a greedy algorithm was used. Note, compatibility 
estimates between design alternatives for system components are not considered. 
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Table 8.6. Design alternatives, estimates and priorities 



Design 
alternatives 


Criteria 


Priority 

n 


Ti 


T 2 


T 3 


T 4 


Bi 


Q 
O 


z 


1 
1 


z 


2 


B 2 


z 




z 


1 
1 




z 


1 


B 3 


A 

4 


1 
1 


1 
1 


1 
1 


2 


E/i 


Q 




z 


1 
1 




z 


2 


u 2 


z 




z 


1 
1 






1 


u 3 


A 

4 


1 
1 


1 
1 


1 
1 


2 


Ri 


1 
1 


Q 




1 
1 


Q 




3 


R 2 


Z 




z 


z 


Q 
O 


3 


R 3 


z 




z 




z 


z 


1 


i?4 


Q 
O 


1 

1 







1 

1 


2 


Vi 


1 
1 




z 


1 
1 


Q 
O 


3 


v 2 


z 


1 
1 


1 
1 


z 


2 


v 3 






1 
1 


z 


1 
1 


1 




4 


1 


2 


1 


2 




2 


1 


2 


1 


1 







2 


1 


2 


2 


Oi 


2 


2 


1 


2 


1 


o 2 


3 


1 


2 


2 


2 




1 


1 


3 


1 


2 




2 


2 


1 


3 


2 


F 2 


3 


1 


2 


2 


1 




1 


1 


1 


1 


1 


D 2 





2 


2 


2 


2 


Gi 


4 


1 


1 


1 


1 


G 2 


4 


2 


3 


2 


3 


G 3 


3 


2 


3 


2 


3 


G 4 


2 


3 


2 


3 


2 




Fig. 8.14. Solution S 2c 
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Table 8.7. Modification operations 



Operations 


Binary 


Cost 


Priorities 




Veil IcLUlC 


a ij 




1. Addition 








1.1. A 1 


Xn 


1 


1 


1.1. None 


Xl2 





3 


1.2. Pi 


X21 


3 


2 


1.2. None 


X 2 2 





3 


1.3. Li 


X31 


1 


1 


1.3. None 


X 3 2 





3 


2. Deletion 








2.1. £a 


X41 


1 


1 


2.2. None 


X42 





2 


3. Replacement 








3.1. B 2 => B 3 


X 5 1 


4 


1 


3.1. None 


X52 





3 


3.2. C/ 2 C/i 


Xei 


3 


2 


3.2. None 


X62 





3 


3.3. Oi => 3 


Xji 


1 


1 


3.3. None 


X72 





3 



8.4. Integrated Security System. Here a numerical example for integrated se- 
curity system is considered from |117) . The initial morphological structure is the 
following (Fig. 8.15): 

0. Integrated security system S = A* B * O. 

1. Closed-circuit television (CCTV) A = J * D: 

1.1. Cameras J: conventional J\ (2), "Tirret" Ja(2), "varifocal" Js(3), and "auto- 
house" J 4 (l). 

1.2. Lighting D: natural -Di(3), natural and guard £>2(2), and natural, guard, 
and alarm D 3 (l). 

2. Access control B = G*U *V: 

2.1. Access to territory G: card Gi(l), radio-pendant 6*2(3), and biometry G 3 (2). 

2.2. Access to building U: card Ui(l), radio-pendant i7 2 (3), and biometry 1/3(2). 

2.3. Access to premises V: card Vi(l), code ^(2), biometry V 3 (2), V4 = 
VikV 2 {2), V 5 = VikV 3 (3), and V 6 = Vx^kV^?,) . 

3. Burglar alarm O — X * Y * Z: 

3.1. Border-line 1 based on some principles X: single physical principle Xi(l), 
two physical principles A 2 (2), and three physical principles A 3 (3). 

3.2. Border-line 2 based on some principles Y: single physical principle li(l), 
two physical principles ^(3), and three physical principles ^3(2). 

3.3. Border-line 3 based on some principles Z: single physical principle Zi(l), 
two physical principles Z 2 (3), and three physical principles Z 3 (2). 

In |117j . two resultant solutions have been obtained (via HMMD): 

5 1 = A 1 *Bi*0 1 = (J 2 *Di)*(Gi*Ui*Vi)*(X 3 *Y 3 *Z 3 ), 

5 2 = A 1 *B 1 *0 1 = ij 2 *Di) + (G 1 *U 1 *Vi)*ix 1 *Y 1 *Z 1 ). 
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,0 = X*Y*Z 

1 = x 3 *y 3 *z 3 

2 = Xi * Yi * Z 1 



X Y 



Z 



G 



S = A*B-kO 

51 = A 1 *B 1 * Oi 

5 2 = A 1 *B 1 * 2 



B = G*U*V 
B 1 = G 1 *U 1 * Vi 



U 



V 



A = J*D 
Ai=J 2 * D 1 



6 



J 



D 



Ji(2) £>i(3) 
J 2 (2) #2(2) 
J 3 (3) D 3 = D x kD 2 {l) 



o o o o 

X 1 {1) Yt{l) Zi(l) Gi(l) E/i(l) Vi(l) 
X 2 (2) y 2 (2) Z 2 (2) G 2 (3) E/a(3) F 2 (l) 
X 3 (3) F 3 (3) Z 3 (3) G 3 (2) (7 3 (2) V 3 (2) 

y 4 = ^i&F 2 (2) J 4 (l) 

14 = Vi&V 3 (3) 

y 6 = vi&F 2 &y 3 (3) 

Fig. 8.15. Basic hierarchy for integrated security system [117] 
In Fig. 8.16 and Fig. 8.17, supersolution and subsolution are depicted. 



J D G U V X Y 



X 3 



ft 

Y 3 



Z 

z 3 



Fig. 8.16. Supersolution 
JDGUyXYZ 

^Cd^(g^(^(v?^^^ 



Fig. 8.17. Subsolution 9 

Here, the obtained subsolution can be used as "system kernel", i.e., K = 0. 
Further, the addition strategy as simplified multiple choice problem is used: 



3 2 

s - t -z2z2 a 

i=i j=i 



max im 

2 



m = l,3;xy e {0,1}. 



A list of addition operations is presented in Table 8.8 |117j . Here cost estimates are 
based on expert judgment, priorities {r^} were computed via Electre-like method 
in [11 7) and transformed to {ci}, b = 7.00. 

A resultant solution is depicted in Fig. 8.18 (xn — 1, xyi — 0, x%\ — 1, a; 22 = 0, 
.t 3 i = 1, .t 32 = 0). Note, compatibility estimates between design alternatives for 
system components X, Y, Z are not considered. 
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Table 8.8. Addition operations 



Improvement 


Binary 


Cost 


Priorities 


operation 


variable 




n 


x x 


Xll 


3 


1 


x 3 


Xl2 


4 


3 




X21 


2 


1 


Y 3 


%22 


3 


3 


Zi 


a?3i 


2 


1 


z 3 


^32 


3 


3 



J D G U y X Y Z 




Fig. 8.18. Resultant aggregated solution 



8.5. Common Educational Course. Here the following initial sets are consid- 
ered (Table 8.9): 

(i) initial set of educational modules for a basic course on combinatorial opti- 
mization A = {1, i, n} = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16}, 

(ii) initial set of educational modules for a course on combinatorial optimization 
for students in "communication systems" A 1 = {1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 12, 13, 14} C 
A, and 

(iii) initial set of educational modules for a course on combinatorial optimization 
for students in "information systems" A 2 = {1, 2, 3, 4, 5, 6, 7, 9, 10, 12, 15} C A. 

The aggregation problem is: 

Design a common course for students in "communication systems" and "infor- 
mation systems" A C A while taking into account weights of educational modules 
({ Wi }, {w j}, {w f} ). 

Here two solving strategies are considered: 

1. Addition strategy (strategy I). 

2. Median-based strategy. 

Let us examine the first case. The following "system kernel" is considered: 
K = A 1 P| A 2 = {1, 2, 3, 4, 5, 7, 9, 10, 12}. The set of elements for addition is: 



B = {i e A\(A 1 f]A 2 )\(w} > 0.73) U (w 2 > 0.73)} = {6,8,11,14,15}. 

Here, three design versions (alternatives) are considered i G B: VI (None), V 2 l 
(compressed version), and V£ (normal version) (j — 1,3). Table 8.10 contains the 
description of the addition elements and corresponding versions. 
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Table 8.9. Illustrative numerical example 



2 


Topic 


A 




A 1 




A 2 


w( A 2 ) 


A 1 II A 2 


A 1 HA 2 
i i 


l 


Optimization 


* 


1.0 




0.6 




0.6 




■k 


2 


Complexity issues 




1.0 


•k 


0.7 


-k 


0.7 


•k 


-k 


3 


Selection/sorting 


-k 


1.0 


-k 


0.6 


•k 


1.0 


-k 


-k 


4 


Knapsack problem 


•k 


1.0 


-k 


0.4 


-k 


1.0 


•k 


-k 


5 


Clustering 


* 


1.0 


* 


1.0 


-k 


1.0 


-k 


■k 


6 


Multiple-choice 
problem 




1.0 




0.0 


-k 


0.76 


* 




7 


Spanning tree 




1.0 


★ 


1.0 


-k 


0.3 


* 


■k 


8 


Routing 




1.0 




1.0 




0.0 






9 


Assignment 




1.0 




0.8 


•k 


0.4 




-k 


10 


Scheduling 




1.0 




0.4 


•k 


0.8 




-k 


11 


TSP 


* 


1.0 




0.75 




0.0 






12 


Covering 




1.0 




0.8 


-k 


0.7 


★ 


■k 


13 


Steiner tree 


* 


1.0 




0.7 




0.0 






14 


Graph coloring 


* 


1.0 




0.73 




0.0 






15 


SAT problem 


* 


1.0 




0.0 


-k 


0.8 






16 


Alignment 




1.0 




0.0 




0.0 







Table 8.10. Addition elements 



K 


Binary 
variable 


Version 

v i 


i 


c Kj 


a Kj 


1 


x n 


V? 


6 


0.0 


0.0 


1 


X\2 


v 2 6 


6 


0.4 


0.5 


1 


£13 


vi 


6 


0.76 


1.0 


2 


X21 


v? 


8 


0.0 


0.0 


2 


X22 


vi 


8 


0.5 


0.5 


2 


X23 


vi 


8 


1.0 


1.0 


3 


X31 




11 


0.0 


0.0 


3 


X32 


v 2 n 


11 


0.45 


0.5 


3 


%33 


v 3 n 


11 


0.75 


1.0 


4 


X41 




14 


0.0 


0.0 


4 


X42 




14 


0.4 


0.5 


4 


X43 




14 


0.73 


1.0 


5 


X51 




15 


0.0 


0.0 


5 


X52 


vi* 


15 


0.4 


0.5 


5 


X53 


vi 5 


15 


0.8 


1.0 



Thus, the addition problem (multiple choice problem) is: 

5 3 

C-KjX K j 

K =i j=i 

53 3 
s-t. ^2 ^2 a Kj x K <b,^2 x *j ^ 1Vk > x kj e {°. !}• 

«=1 j=l j=l 
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Here c K 3 = max(w^, w 2 ), estimate of a K j is based on expert judgment (Table 8.10), 
and b = 3.50. 

A resultant solution based on strategy I is depicted in Fig. 8.19 (xu = 0, 

X12 = 0, Xi 3 = 1, X 2 1 = 0, X 2 2 = 0, X 2 Z = 1, £31 = 1, X 32 = 0, X33 = 0, X41 = 0, 

X42 = 1, X43 = 0, .T51 = 0, X52 = 0, .T53 = 1). Evidently, normal version V3 is 
used for elements of "system kernel" K. The solving process was based on a greedy 
algorithm. Here compatibility estimates between design alternatives for system 
components were not considered. 



6 7 



10 11 12 13 14 15 16 



c*r\ ^ c5\ <5\ ^ & cZ\ ^\ c*^ c*^ c*^ ^ 



v 



vi 



v 3 12 



V, 



15 



Fig. 8.19. Aggregated common course (addition strategy) 

In the second case, Pareto-efficient median-solutions {A } are searched for through 
the following two-component vector criterion p = (p 1 (A°, A 1 ), p 2 (A°, A 2 )) (by weights 
{w}} and {w 2 }): 

min^^U 1 ) = 1 - ^ A °^ Al \ W l minp 2 (A^A 2 ) = 1 - fi^H^ . 

l^ie(A°\JA 1 ) W i l^ie(A a [\A 2 ) W i 

Generally, the problem of searching for the median belongs to class of NP-hard 
problems. Let us consider an approximation heuristic. Initial elements for the 
median are the following (i.e., deletion of element 16): 



{1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15}. 



Thus, cardinality of search space is: = 2 15 . 

First, let us assume that elements {1,2,5} will be included into each median- 
solution (to decrease problem dimension) and, as a result, {A } = 2 12 . 

Second, elements from set {3, 4, 9, 10, 11, 12} will be included into each median 
solution as well and, as a result, |{^4°}| = 2 6 = 64. 

Finally, only six submcdian-subsolutions are selected for examination: 

A 01 = {7, 8, 13, 14, 15}, A 02 = {6, 8, 13, 14, 15}, A 03 = {6, 7, 13, 14, 15}, 

A 04 = {6, 7, 8, 14, 15}, A 05 = {6, 7, 8, 13, 15}, A 06 = {6, 7, 8, 13, 14}. 

Table 92 contains the subsolutions (k = 1,6), corresponding values of vector 
~p(A° K ) = (p lK (A° K , A 1 ), p 2k (A 0k , A 2 )), and information on inclusion into the layer 
of Pareto-efficient solutions. Here parts of A 1 and A 2 are considered, which corre- 
spond to {6, 7, 8, 13, 14, 15}: {7, 8, 13, 14} (for A 1 ) and {6, 7, 13, 15} (for A 2 ). 

Thus, the following preference relation is obtained over the six subsolutions 
above: A 01 h A 06 and A oi h A 05 h A 03 h A 02 . 

Finally, subsolutions A 01 and A 04 are Pareto-efficient ones (Table 8.11, Fig. 
8.20). 



AGGREGATION OF COMPOSITE SOLUTIONS: STRATEGIES, MODELS, EXAMPLES 47 



Table 8.11. Candidate-subsolutions 



K 


Candidate subsolution 
A° K 


p 1 (A° K ,A 1 ) 


p 2 (A° K ,A 2 ) 


Inclusion into 
Pareto-layer 


1 


A 01 = {7,8,13,14,15} 


0.00 


0.40 




2 


A 02 = {6,8,13,14,15} 


0.30 


0.05 




3 


A 03 = {6,7,13,14,15} 


0.30 


0.00 




4 


A 04 = {6,7,8,14,15} 


0.21 


0.00 




5 


A 05 = {6,7,8,13,15} 


0.22 


0.00 




6 


A oe = {6,7,8,13,14} 


0.00 


0.43 





P 2 (A°,A 2 ) 
1 



G) 



A 01 

The point of ideal 
solution(s) 



(0,0)4 

Fig 



A 



04 



p\A°,A^ 



.20. Pareto-efficicnt solutions 



The resultant (extended) Pareto-efhcient solutions A 01 and A 04 are presented in 
Fig. 8.21 and in Fig. 8.22 (normal versions of course components are used). 



1 
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vi 
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vi 
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15 



Fig. 8.21. Median-based aggregated common course A 01 
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Fig. 8.22. Median-based aggregated common course A 04 



8.6. Plan of Students Art Activity. The considered numerical example is a 
small part of the example of students plan from [101 . The initial morphological 
structure is the following (Fig. 8.23): 

1. Plan of students art activity S = I * J -kU: 

1.1 dance I: I\ (none), I 2 (ball dance), 23 (ensemble); 

1.2 music J: J\ (none), J 2 (classic), J 3 (jazz), J4 (singing); 
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1.3 theatre U: U\ (none), U 2 (actor), U3 (producer), U4 (technical worker), U5 
(author) . 

The following criteria were used for assessment of the design alternatives: cost 
Ci, opportunity to meet new friends C 2 , opportunity to meet boy friend or girl 
friend C3, accordance to personal inclinations C4, usefulness for future career C5, 
usefulness for health Cq, and usefulness for future life C7. The resultant priorities 
of design alternatives (via Electre-like technique) are shown in parentheses (Fig. 
60). Thus, three resultant solutions are the following (via HMMD): 

Si = h * h * U2, S 2 = h * ^3 * U 2 , and S3 = I 3 * J 2 * U 5 . 

In Fig. 8.24 and Fig. 8.25, supersolution and subsolution are depicted. Note, 
the subsolution contains only one element. Thus, the design of "system kernel" 
is based on the special method with a = 0.6 (Fig. 8.26). Finally, the obtained 
"system kernel" K can be considered as the resultant solution. 



S = I*J*U 

51 = h * J3 * U 2 

5 2 = h * ^3 * U 2 

5 3 = h * J 2 * U 5 



I J 



o 
h (3) 

Hi) 



u 



■A(4) 
J a (2) 
^3(2) 
J 4 (3) 



U 2 (l) 
U 3 (2) 

U 5 (l) 



J u 



■h 

■h 



U 2 

u 5 



Fig. 8.23. Basic hierarchy Fig. 8.24. Supersolution 



r J U 

13 



I J u 



Fig. 8.25. Subsolution 6 Fig. 8.26. "System kernel" K 



8.7. Combinatorial Investment. The considered numerical example was pre- 
sented in |101j . The initial morphological structure is the following (Fig. 8.27): 



S = A* B * L 
S"i = A 4 * B 3 * Li 
S 2 =A 2 *B 5 * Li 
S 3 =A 2 *B 3 * Li 

1S4 = A 2 -k B5 -k L4 



6 A B 

Ax (3) 
A 2 (3) 

Ms) 

Ai(2) 



L 



Br (3) 
B 2 (4) 
B 3 (2) 
B 4 (3) 
B 6 {1) 



Lx(3) 
£ 2 (3) 
£ 3 (2) 
L 4 (l) 



A B L 



A 2 

A4 



B,, 
B, 



Fig. 8.27. Basic hierarchy Fig. 8.28. Supersolution 
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1. composite portfolio S = A* B ★ L: 

1.1 short-time investment A: A\ (state bonds), A2 (bank deposit), A3 (specula- 
tion on the stock exchange), A4 (oil shares); 

1.2 middle-time investment B: B\ (state bonds), B2 (bank deposit), B 3 (im- 
movables), B4 (jewelry), B 5 (shares in biotechnology); 

1.3 long-time investment L: L± (state bonds), L2 (bank deposit), L 3 (antique), 
L4 (shares of airspace companies). 

The following criteria were used for assessment of the design alternatives: pos- 
sible profit Ti, risk T2, prestige T3, possibility for continuation T4, possibility to 
establish a new company T5, obtaining a new experience Tg, possibility to organize 
a new market T7, possibility to obtain "name" Ts, and connection with previous 
activity T9. 

The resultant priorities of design alternatives (via Electre-like technique) are 
shown in parentheses (Fig. 8.27). Thus, four resultant solutions are the following 
(via HMMD): 

Si = A4, *B 3 ★la, S 2 = A 2 *B 5 ★ii, 5*3 = A 2 ★ B 3 ★ L 4 , and S4 = A 2 ★ B 5 ★ L 4 . 

In Fig. 8.28 and Fig. 8.29, supersolution and subsolution are depicted. Note, the 
subsolution does not contain elements (empty). Thus, the design of "system kernel" 
is based on the special method with a = 0.6 and selection of the best elements for 
system parts: B$ for components B and L4 for components L (Fig. 8.30). Finally, 
the obtained "system kernel" K can be considered as the resultant solution. 



A B L 



A B 1 



Fig. 8.29. Subsolution 6 Fig. 8.30. "System kernel" K 



8.8. Configuration of Applied Web-based Information System. Here ag- 
gregation of solution for applied Web-based information system is considered |114j . 
The tree-like model of the system is depicted in Fig. 8.31. 

%S = A*B = {J*E)*(W*D*0) 



Hardware 
A = J-kE 



J ^1 

Server Server for 
for DBs applica- 



Software 

B = W*D*0 



■h 
■h 
■h 



tions 
Ei 
E% 
E3 
Ei 



W 
Web 
server 
Wi 
W 2 
W 3 
W 4 
W 5 



DBMS 
Di 
D 2 
D 3 



"lo 

Operation 
system 

Oi 

O2 

Ol 

o 5 



Fig. 8.31. Structure of Web-based system [1 14] 



DAs for system components are the following jl 14] : 

(1) server for DBs J: PC (Ji), Supermicro (J 2 ), and Sun (J 3 ); 

(2) server for applications E: on server of DBs (£1), Sun (E 2 ), Supermicro (£3) 
and PC {Ei)- 



50 



MARK SH. LEVIN 



(3) Web-server W: Apache HTTP-server (Wi), Microsoft IIS (W 2 ), Bea We- 
blogic (W3), Web Sphere (W4), and Weblogic cluster (W5); 

(4) DBMS D: Oracle (£>i), Microsoft SQL (D 2 ), and designed SQL (D 3 ); and 

(5) operation system O: Windows 2000 server (Oi), Windows 2003 (0 2 ), Solaris 
(O3), FreeBSD (0 4 ), and RHEL AS (0 5 ). 

In |114j . the design problem was solved for three basic applied situations: (a) 
communication provider (Fig. 8.32), (b) corporate application (Fig. 8.33), and (c) 
academic application (Fig. 8.34). 



• S=A*B 
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E 4 {3) W 5 {3) 6 (1) 
Fig. 8.32. Communication provider [1 14] 
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Fig. 8.33. Corporate application |114) 
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I S = A-kB 

Si =A\*B\ = (J 3 * E 2 ) * (Wi * £ 2 * 3 ) 
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Fig. 8.34. Academic application [114) 



Now let us consider an aggregation problem: 

To obtained an aggregated solution while taking into account all three application 
situations. 

Here two simplified heuristic solving strategies are examined: 

(i) aggregation of composite solutions for three considered cases; 

(ii) aggregation of information at the first stage of the solving process (aggre- 
gation of rankings for DAs of system components J, E, W, D, O) and usage of 
HMMD to solve the design problem. 

Let us describe the first strategy. Here the solving process is based on median/kernel- 
based strategy (Fig. 8.35). The resultant solution is (an approximate median/kernel): 
S = (J 2 *E 2 )*(W 1 *D 3 *0 2 ). 



Solutions for provider 

= J 2 *E 2 *W 1 *D 3 *0 3 
J 2 * E 2 * W 2 * D 2 * 2 
J 2 -kE 2 -kWi*D 2 -kO, 



Aggregation process 




S a " = J 2 *E 2 *W 1 *D 3 *0 2 



Solutions for corporate application 

S? = Ji * E x * Wi * D 3 * 5 
Si = Ji *E 1 -kW 2 -kD 3 -k0 2 
Si = J 2 *E 2 *W 1 *D 3 *0 5 
01 = J 2 *E 2 *W 2 -kD s -k0 2 

Solutions for academic application 
J 3 * E 2 * Wi * D 2 * 3 



Fig. 8.35. Aggregation of solutions 



The second solving strategy is based on preliminary aggregation of rankings (Fig. 
8.36, Fig. 8.37, Fig. 8.38, Fig. 8.39, Fig. 8.40) and a new design. 
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Fig. 8.38. Aggregation for W Fig. 8.39. Aggregation for D 
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Fig. 8.40. Aggregation for O 
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After the usage of HMMD for the aggregated rankings of DAs (for system com- 
ponents J, E, W, D, O), the following four resultant composite DAs are obtained 
(Fig. 8.41, compatibility estimates are contained in |114j ): 

(a) Si" =A 1 *B 1 = (J 2 *E 2 )* (Wi * D 2 * 5 ), 

(b) S* 99 = A X *B 2 = (J 2 * E 2 ) * (Wi *D 3 * 8 ), 

(c) S%" = A 2 *B 1 = ( J 3 * E 2 ) * (Wi * D 2 * 5 ), and 

(d) S 4 a9ff = A 2 * B 2 = ( J 3 * £ 2 ) * (Wi * £> 3 * 5 ). 



5 = ,4*5 

5"" =A 1 + B 1 = (J 2 * E 2 ) * (Wi *D 2 * B ) 
Ai * S 2 = ( J 2 * £ 2 ) * (Wi * £> 3 * 5 ) 
A 2 * Si = ( J 3 * £2) * (Wi * D 2 * 5 ) 

^2*S 2 = (J 3 *£ 2 )*(Wi*L> 3 *0 5 ) 

Hardware Software 

v4 = J*£ B = W*D*0 
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J 3 (2) £ 2 (1) W 3 (3) 3 (1) 

£ 3 (2) W 4 (3) 4 (2) 

£ 4 (3) W 5 (3) 5 (1) 

Fig. 8.41. Aggregation as design (2nd strategy) 



8.9. Modular Educational Course on Design. Here aggregation of three ed- 
ucational courses (morphological structures) is examined: (1) course on systems 
engineering (structure A 1 ) ( |101j . [102j ) (Fig. 8.42); (2) course on information en- 
gineering (structure A 2 ) ([99], [101] ) (Fig. 8.43); and (3) course on hierarchical 
design (structure A 3 ) ([M], [106], p33]) (Fig. 8.44). The following general types 
DAs (i.e., of the corresponding educational module) are examined for each leaf 
node of the presented hierarchical models (Fig. 8.42, Fig. 8.43, and Fig. 8.44): 
"absence" of the educational module Xq, compressed information on the educa- 
tional module Xi, teaching at a medium level X 2 , serious teaching X 3 , and serious 
teaching with a special student research work/project X4. Compatibility estimates 
for three examined morphological structures above are presented in Tables (expert 
judgment): (i) structure A 1 : Table 8.12, Table 8.13; (ii) structure A 2 : Table 8.14, 
Table 8.15, Table 8.16; and (iii) structure A 3 : Table 8.17, Table 8.18, Table 8.19. 

After usage of HMMD, the following composite DAs are obtained: 

I. For structure A 1 (Fig. 8.42): 

D 1 = L 3 *G 2 *C 2 *M 2 , JV(Di) = (2; 3, 0,1), 

D 2 = L i *G 2 *C 1 *M 2 , N{D 2 ) = (1;2,2,0); 

Qi = A 3 *B 2 , N(Qi) = (3; 2, 0, 0), Q 2 = A 4 * B 4 , N(Q 2 ) = (3; 2, 0, 0); 
S{ =D 1 *Q 1 = (L 3 *G 2 *C 2 * M 2 ) * (A 3 * B 3 ), 
S\ =D 1 *Q 2 ^ (L a * G 2 * C 2 * M 2 ) * (A 4 * B 4 ), 
Si =D 2 irQ 1 ^ (L 4 * G 2 * Ci * M 2 ) * (A 3 * B 3 ), 
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Si =D 2 *Q 2 = (Li * G 2 * & * M 2 ) * (Ai * B 4 ). 

II. For structure A 2 (Fig. 8.43): 

h = L 3 * G 2 * G 2 * £ 2 , JV(Ji) - (3; 2, 1, 0); 

Oi = # 3 *W 2 *M 2 , AT(Oi) = (2;3,0,0),O 2 = tf 2 *W 2 *M 2 , AT(0 2 ) = (3; 2, 1,0); 
Qi = K 3 * J 3 * ^2 * i? 3 , iV(Qi) = (3; 4, 0, 0); 

Si =I 1 *0 1 *Q 1 = (G 2 * G 2 * £ 2 ) * (H 3 *W 2 * M 2 ) * (V 3 * J 3 * F 2 * R 3 ), 
S$ = h*0 2 *Qi = (G 2 *C 2 * E 2 ) * (H 2 *W 2 * M 2 ) * (V 3 * J 3 * Y 2 * i? 3 ). 

III. For structure A 3 (Fig. 8.44): 
h = G 2 *C 2 , N(h) = (3; 2, 0,0); 

01 = H2*W2*M 2 *U lt JV(Oi) = (3; 3, 1,0), 

2 = i? 3 * VF 2 * M 2 * C/i, JV(0 2 ) = (2; 4, 0, 0); 
Qi = F 4 * Ji * y 2 * Z 3 , JV(Qi) = (3; 3, 1,0); 

Sf =I 1 *Oi*Q 1 = (G 2 * C 2 ) * (JT 2 * W 2 * M 2 * C/i) * (F 4 * Ji * F 2 * Z 3 ), 

Sf = A * 2 ★ Qi = (G 2 * G 2 ) * (H 3 * VF 2 * M 2 * C/i) * (F 4 * Ji * F 2 * Z 3 ). 

Evidently, it is possible to aggregate the obtained composite solutions (Fig. 8.45). 
On the other hand, let us consider the following extended aggregation strategy IV. 
General structure (A) consists of the following parts: (i) tree-like system model T, 
(ii) set of leaf nodes P, (hi) sets of DAs for each leaf node D, (iv) DAs rankings 
(i.e., ordinal priorities) R, and (v) compatibility estimates between DAs I. Thus, a 
vector proximity for two structures A Q , A^ can be examined as follows: 

p (A", A") = (p*(T a ,T /3 ), p*(P«,P^), p*(D a ,D^), /j t (R Q ,R /3 ), P \l a ^)). 

As a result, A a " has to be searched for as Pareto-efficient solution(s) by the 
vectors ~p(A a " , A 1 ) Mi e {i} where index i corresponds to an initial solution. This 
problem is very complicated. Let us consider a simplified solving framework: 

Phase 1. Aggregation of basic initial data for initial structures: 

1.1. aggregation of morphological structures including the following: 
(1.1.1.) aggregation of sets of leaf nodes, 

(1.1.2.) aggregation of sets of DAs for each leaf node, 

(1.1.3.) aggregation of DAs rankings, and 

(1.1.4.) aggregation of compatibility estimates for DAs sets; 

1.2. aggregation tree- like structures. 
Phase 2. New hierarchical design. 

Thus, the following stages are considered for examined three structures (A 1 , A 2 , 
A 3 ): 

Stage 1. Aggregation of leaf node sets for the initial structures (Fig. 8.46); 
Stage 2. Aggregation of morphological structure: 

2. 1 . aggregation of sets of DAs for each leaf node (by each leaf node (while taking 
into account addition of DAs which correspond to "absence", i.e., aggregation of 
rankings with extension of DAs sets) (Fig. 8.47, Fig. 8.48, Fig. 8.49, Fig. 8.50, 
Fig. 8.51, Fig. 8.52, Fig. 8.53, Fig. 8.54, Fig. 8.55, Fig. 8.56, Fig. 8.57, Fig. 8.58, 
Fig. 8.59, Fig. 8.60, Fig. 8.61, Fig. 8.62); 

2.2. aggregation of interconnection (compatibility estimates) for DAs sets (Ta- 
bles 8.20, 8.21, 8.22, 8.23; selection of a maximal value or expert judgment); 

Stage 3. Building of an aggregation tree-like structure (Fig. 8.63); and 
Stage 4- Hierarchical design of an integrated course (as an aggregated solution) 
(Fig. 8.63). 

Here the aggregated superstructure A (Fig. 8.63) has been obtained via expert 
judgment. 
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S 1 = D * Q 

Sj =£>i*Qi = (L 3 *G 2 *C 2 *M 2 )*(A 3 *B 3 ) 
D 1 *Q 2 = (L 3 * G 2 * G 2 * M 2 ) * (At * B 4 ) 
(L 4 * G 2 * Ci * M 2 ) * (A3 * B3) 
(Li * G 2 * Ci * M 2 ) * (A* * Bi) 
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Fig. 8.42. Course on systems engineering (structure A 1 ) 
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Fig. 8.43. Course on information system (structure A 2 ) 
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S 3 = I * O * Q 
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Fig. 8.44. Course on hierarchical design (structure A 3 ) 
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Fig. 8.45. Aggregation of solutions 
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Table 8.12. Compatibility 
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Table 8.13. Compatibility 
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Table 8.15. Compatibility 
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Table 8.16. Compatibility 
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Table 8.17. Compatibility 
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Table 8.18. Compatibility 
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A 2 




Fig. 8.46. Aggregation of leaf node sets 
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Fig. 8.47. Aggregation for L 






Fig. 8.48. Aggregation for G 
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Fig. 8.49. Aggregation for C 
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Fig. 8.50. Aggregation for E 
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Fig. 8.52. Aggregation for W 
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Fig. 8.55. Aggregation for V Fig. 8.56. Aggregation for F 
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Fig. 8.61. Aggregation for A 
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_S = I-kO-kX 

S 1 = I 1 *0i* Xi 

s 2 = h * Oi * Xi 

s 3 = / 3 *Oi*X 1 
Methodology 
/ = LkGkCkE 
I\ = L 3 ★ G 2 * G 2 * E 2 

12 = L 3 -k G 2 -k C 3 -k E2 

1 3 = L4 * G3 * C3 * E 2 
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Fig. 8.63. Superstructure for aggregated course (structure A) 
Table 8.20. Compatibility 
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Table 8.21. Compatibility 
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Table 8.22. Compatibility 
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Table 8.23. Compatibility 
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After usage of HMMD, the following composite DAs are obtained for the resul- 
tant aggregated structure (A) (Fig. 8.63): 

Xi = V 3 *F 4 *J 3 *R 3 *Y 2 *Z 3 *A 3 *B 3 , N(Xi) = (3; 8, 0,0), 
O x = H 3 * W 2 * M 2 * U u JV(Oi) = (2; 4, 0, 0); 
h = L 3 *G 2 *C 2 *E 2 , N(Oi) = (2; 4, 0,0), 
h = L 3 *G 2 *C 3 *E 2 , N(Oi) = (2; 4, 0,0), 
h = U * G 3 * G 3 * E 2 , N(Oi) = (3; 2, 2, 0); 
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S 1 =I 1 *0 1 *X 1 = (L 3 * G 2 * C 2 * E 2 ) * (H 3 *W 2 * M 2 * U x ) * (V 3 * F A * J 3 * 
R 3 -kY 2 -kZ 3 -kA 3 -kB 3 ), 

S 2 = I 2 *Oi*Xi = {L 3 *G 2 *C 3 *E 2 )*(H 3 *W 2 *M 2 *U 1 )*(V 3 *F 4 *J 3 * 
R 3 *Y 2 * Z 3 *A 3 * B 3 ), 

S 2 = I 3 * Oi * Xi = (£3 * G 3 * C 3 * £2) * (#3 * ^2 * M 2 * C7i) * (V3 * F A * J 3 * 
^3*^2*^3*^3*^3). 



8.10. Configuration of Car in Electronic Shopping. Recently, many products 
have a complex configuration and a buyer can often generate a product configura- 
tion that is more useful for him/her. An approach to electronic shopping based on 
product configuration (morphological approach) has been suggested in [109 . Here 
let us consider a multi-choice scheme for selection of structured product /system in 
electronic shopping with a resultant aggregation (Fig. 8.64). 
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Fig. 8.64. Multi-choice scheme of electronic shopping 

Now the following example is described (at a conceptual level). An initial 
morphological structure A of a car is the following (Fig. 8.65) (in real applica- 
tion, this structure can be considered as a result of processing the selected prod- 
ucts/solutions): 

0. Car S = A-kB-kC. 

1. Main part A = E * D: 

1.1. Engine E: diesel E\, gasoline E 2 , electric E 3 , hydrogenous E4, and hybrid 
synergy drive HSD £5; 

1.2. Body D: sedan D\, universal D 2 , jeep D 3 , pickup D 4 , and sport D 5 . 

2. Mechanical part B = X *Y * P * Z: 

2.1. gear box X: automate X±, manual X 2 ; 

2.2. suspension Y: pneumatic Y\, hydraulic Y 2 , and pneumohydraulic Y 3 ; 

2.3. drive P: front-wheel drive P\, rear-drive P 2 , all- wheel-drive P 3 . 

3. Safety part C = O * G: 

3.1. O: "absence" Oo, electronic Oi; 



64 



MARK SH. LEVIN 



3.2. Safety subsystem G: "absence" Go, passive Gi, active G 2 . 
• A = A*B*C 
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A = E*D 
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Fig. 8.65. General structure of car A 

The following initial solutions (obtained by users) are considered: 
Si =E 1 *D 1 *X 1 *Y 1 *P 1 *Oi*Gi, S2 1 = E 5 *D 1 *X 1 *Y 1 *P 1 *0 1 * G 2 , 
Sf = E 2 *D 1 *X 2 *Y 1 *Pi*O *Gi, Sf = E2*D 3 *X 1 *Y 2 'kP3-kOi*G , and 
S| = £ 2 * £5 * ATi * Y 3 * Pi * Oi * Gi. 

An aggregated solution can be considered as simplified aggregation of five sets 
above: 

S agg =E 2 *D 1 *X 1 *Y 1 *P 1 *0 1 *G 1 . 



9. Conclusion 

In the article, a systemic view point to aggregation of modular solutions is firstly 
presented. The aggregation strategies have been suggested and considered: (i) ex- 
tension strategy (addition of elements to a "system kernel"), (ii) compression strat- 
egy (deletion of elements of a superstructure), (iii) combined extension/compression 
strategy (addition, deletion, and replacement of elements in "system kernel" ) , and 
(iv) strategy of a new design (while taking into account new elements). The exam- 
ined strategies have not been previously discussed in the literature. Our material 
corresponds to the first step in the investigated domain (aggregation of modular 
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solutions). The considered aggregation approaches are useful for engineering do- 
mains, management, computer science and information technology. 

It is necessary to point out, close aggregation problems have been intensive stud- 
ied and used in several domains: (a) decision making (aggregation of rankings to 
obtain a consensus, aggregation of preferences) (e.g., [35], [36], [166] ); (b) integra- 
tion of organizational structures in organization science (e.g., [10], [38], [39], [41] . 
[58)): and (c) integration of information (database schema integration, integration 
of knowledge base structures, integration of catalogs, merging and integration of 
ontologies) (e.g., [JJ, [H], [32], [123], [T3T], [131], [159] V 

In the article, many illustrative applied examples are presented. It is neces- 
sary to point out fundamentals of the described approaches correspond to practice 
(e.g., engineering, management). Thus, the suggested composite solving schemes 
and their local parts have to be based on a combination of mathematical com- 
binatorial models and expert judgment of domain experts. Our approach for 
modeling the complex modular solutions/systems is based on the following four 
layers: (i) system hierarchy as tree-like structure, (ii) system components (leaf 
nodes of the tree- like mode above), (iii) sets of design alternatives for each sys- 
tem components, (iv) compatibility between design alternatives of different sys- 
tem components. Thus, models, problems, and algorithmic schemes are examined 
for the above-mentioned structures: proximities/metrics, substructures (including 
median/consensus, etc.), aggregation strategies/frameworks. The considered solv- 
ing strategies are based on combinatorial problems (multicriteria ranking/selection, 
knapsack problem, multiple choice problem, morphological combinatorial synthesis, 
building median/consensus for structures). 

The suggested aggregation problems may be useful in the following situations: 
(a) design processes based on uncertainty (e.g., fuzzy sets, stochastic models) can 
lead to generating a set (as a grid) of design solutions (here the problem under 
uncertainty can be approximated by a set of deterministic problems and, further, 
their solutions can be aggregated); (b) scenario-based methods can lead to a set of 
design solutions which can be aggregated. Note, complicated methods and models 
are required for dynamical modeling of complex systems, for example: (a) various 
kinds of Petri nets, (b) dynamical graphs. Our material does not involve this type 
of studies. 

In the future, the following research directions can be considered: 

1. usage of the suggested aggregation approaches for various kinds of system con- 
figuration solutions (e.g., set of selected elements, assignment/allocation solutions, 
network-like solutions |110| ): 

2. usage of other types of metrics/proximities for structures; 

3. additional theoretical and applied studies of formal approaches to aggregation 
problems based on aggregation of system tree-like models and/or general system 
structures (i.e., multicriteria searching for median-like tree(s) and/or median-like 
structure(s) as Pareto-efhcient structured solution(s)); 

4. examination of more complicated design problems for several resultant aggre- 
gated solutions (e.g., as "product line"); 

5. examination of various models (for example, building of median/consensus), 
including multicriteria problem statements and usage of special quality domains/ 
spaces, e.g., lattice-like quality domain/space ( [101] . [103] . [106] ): 
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6. examination and usage of other structures for system modeling, e.g., pyramids 

m, my, 

7. taking into account uncertainty; 

8. usage of artificial intelligence methods in problem solving processes; 

9. investigation of new applications; 

10. usage of the suggested aggregation problems as auxiliary underlaying sub- 
problems in information integration/fusion (for example, for approximation of alignment- 
like problems in bioinformatics, image processing; multi-source fusion of structured 
information) ; 

11. special examination of aggregation problems for aggregation of solutions in 
combinatorial optimization problems (e.g., routing, scheduling, traveling salesman 
problem, assignment / allocation, graph coloring, covering) ; 

12. usage of the suggested aggregation approaches in electronic shopping (and 
recommender systems) as the following two-stage framework: (i) an user is selecting 
a preliminary set of modular products (by his/her choice) from product catalogs, 
(ii) a computer system (i.e., special design support service) is constructing an ag- 
gregated product (selection of the best product from the preliminary product set 
is a simplified case); 

13. it may be very interesting to apply the suggested aggregation approaches to 
drug design (generally, to combinatorial chemistry); and 

14. usage of the suggested aggregation approaches in education (engineering, 
management, computer science and information technology). 
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